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Objectives 


The  main  objective  of  this  3-year  project  is  to  grow  high  quality  Bib  single  crystals  (>  1 
cm3  in  volume)  via  a  high  temperature  modified  Bridgman  crystal  growth  technique.  We 
will  then  test  and  analyze  the  crystals  using  a  variety  of  techniques  to  determine  crystal 
quality  (lattice  parameters,  lattice  strain,  impurities,  etc.)  and  potential  methods  to 
improve  Bib  crystals.  Finally,  test  structures  will  be  designed  and  their  performance  will 
be  assessed  using  a  variety  of  small,  calibrated  gamma-ray  sources  available  at  the 
University  of  Florida.  We  will  measure  radiation  sensitivity,  charge  carrier  mobility  and 
lifetimes,  and  noise  characteristics  of  the  test  structures  (basic  material  properties  for 
Bib).  While  the  main  objectives  of  the  project  have  not  changed,  more  emphasis  is 
being  placed  on  Bib  crystal  growth  principles  and  materials  analyses.  This  increase  in 
materials  properties  assessment  is  based  on  suggestions  from  the  DTRA  review  team 
during  the  kickoff  meeting  in  May  2007. 


1 .0  Crystal  Growth  Techniques  of  Bii3 

While  the  PVT  method  can  produce  high  quality  but  small  volume  crystals,  the  vertical 
Bridgman  technique  can  grow  crystals  with  larger  volume  but  the  crystal  quality  is  not  as 
good  as  those  grown  from  PVT.1'4  We  primarily  discuss  the  growth  of  Bil3  single 
crystals  using  the  modified  vertical  Bridgman  (MVB)  technique.  The  quality  of  the 
crystals  grown  under  different  conditions  was  evaluated  and  compared  to  that  of 
crystals  grown  by  PVT  method. 

1.1  Vertical  Bridgman  Technique 

As  the  name  says,  this  technique  is  named  after  P.  W.  Bridgman  who  is  the  first  one 
using  this  method  to  grow  a  series  of  metal  single  crystals.5  The  principle  of  this 
technique  involves  lowering  a  crucible  through  a  furnace  so  that  the  freezing  starts  at 
the  lowest  point  in  the  crucible  and  the  solidification  interface  slowly  moves  up  along  the 
crucible.  The  movement  rates  for  such  processes  range  from  around  0.1  to  200  mm/h.6 
Atypical  vertical  Bridgman  system  is  shown  in  Figure  1.  The  movement  of  the  crucible 
is  controlled  by  a  dropping  motor.  A  longitudinal  temperature  profile  is  established  at  the 
center  of  the  furnace  with  specific  temperature  gradient  near  the  melting  point  of  the 
material.  The  hole  in  the  lid  should  be  small  and  the  lid  should  fit  well  with  the  furnace 
body  to  prevent  thermal  disturbance. 
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Figure  1 .  A  typical  Bridgman  growth  furnace  setup,  (a)  Schematic  diagram  for  the  furnace,  (b) 
Longitudinal  temperature  profile  of  the  furnace. 

The  crystal  growth  with  vertical  Bridgman  method  can  use  a  seed  or  no  seed.  A  non- 
seeded  Bridgman  method  relies  on  random  nucleation  to  produce  a  single  crystal  which 
then  propagates  to  form  one  large  crystal  through  the  whole  frozen  charge.  The 
formation  of  just  one  nucleus  requires  a  small  nucleation  volume.  Thus  the  crucibles  are 
normally  shaped  with  sharp  tips  to  meet  the  nucleation  criteria,  as  illustrated  in  Figure  3- 
1(a).  Fora  seeded  Bridgman  growth,  the  crucible  shape  is  not  that  important.  Normally 
a  flattened  bottom  is  made  to  accommodate  the  seed  material.  The  most  obvious 
requirement  for  a  crucible  is  that  it  should  not  react  with  the  crystals  within  the 
temperature  range  of  crystal  growth.  It  also  should  have  a  smaller  coefficient  of  thermal 
expansion  than  the  crystal  to  avoid  the  induced  strain  and  stress  during  cooling. 

For  vertical  Bridgman  crystal  growth,  the  quality  of  the  single  crystal  is  largely 
determined  by  the  control  of  heat  flow  during  crystallization.  The  thermal  conditions  of  a 
practical  crystal  growth  system  are  rather  complicated,  thus  analytical  solutions  are  not 
possible  without  simplifications  and  assumptions.7  In  the  following  discussion,  it  is 
assumed  that  the  melt  is  a  single  component  material,  the  density  of  the  melt  and  its 
crystal  are  the  same,  temperature  gradients  within  the  crystal  transverse  to  the  growth 
direction  are  negligible,  and  the  solid-liquid  interface  is  planar.  During  crystal  growth,  a 
heat  balance  is  set  up  by  the  conservation  of  heat  flux  at  the  interface,  which  can  be 
written  as8 

PSLR  =  KsGs  -KlGl  (1) 

where  Ks  and  KL  are  the  thermal  conductivity  of  solid  and  liquid  (cal/(cm-°C  s)),  Gs  and 
Gl  are  the  temperature  gradient  in  solid  and  liquid  at  the  solid-liquid  interface  (°C/cm), 
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ps  is  the  density  of  solid  (g/cm3),  L  is  the  latent  heat  of  melt  (cal/g),  and  R  is  the  growth 
rate  (cm/s).  The  left  side  in  Equation  (1)  indicates  the  heat  generated  at  the  interface 
due  to  solidification.  The  first  term  on  the  right  side  is  the  heat  flux  going  out  of  the  solid- 
liquid  interface.  The  second  term  on  the  right  side  refers  to  the  heat  flux  coming  into  the 
solid-liquid  interface.  Equation  (1)  indicates  that  the  growth  rate  R  can  be  controlled  by 
independently  adjusting  Gs  and  GL.  This  is  important  as  growing  good  single  crystals 
requires  that  the  temperature  gradient  be  high  and  the  growth  rate  remain  low.  The 
temperature  gradient  in  liquid  (Gl)  is  directly  proportional  to  the  heat  flux  in  the  liquid  at 
the  solid-liquid  interface.  The  temperature  gradient  in  solid  (Gs)  is  normally  determined 
empirically  or  through  heat  flow  calculations.  In  a  cylindrical  system,  the  heat  conduction 
in  solid  can  be  described  as6 
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where  cs  is  the  specific  heat  of  solid  (cal/(g°C)),  and  z  is  the  coordinate  parallel  to  the 
growth  direction.  Physically,  the  term  on  the  left  represents  the  net  rate  of  heat 
conduction  away  from  the  solid  per  unit  volume.  The  first  term  on  the  right  side 
represents  the  decrease  rate  of  internal  energy  per  unit  volume.  The  second  term  on 
the  right  side  is  the  rate  of  energy  released  per  unit  volume  due  to  interface  moving. 
Some  simplifications  need  to  be  made  to  solve  equation  (2).  For  example,  in  the  case  of 
floating  zone  (crucibleless)  growth,  assuming  the  crystal  has  a  cylindrical  shape, 
temperature  gradients  along  the  transverse  direction  are  negligible,  and  heat  transfer 
from  the  crystal  to  surrounding  is  by  convection.  At  steady  state,  the  temperature 
gradient  in  solid  at  the  solid-liquid  interface  is  given  by8 
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where  T0  is  the  ambient  temperature  (°C),  Tm  is  the  melting  point  of  the  crystal,  as  is  the 
thermal  diffusivity  of  the  solid  crystal  (Ks/pscs,  cm2/s),  h  is  the  heat  transfer  coefficient 
from  solid  to  ambient  (cal/(cm2oCs)),  and  r  is  the  radius  of  crystal  (cm).  Since  the 
growth  rate  R  is  usually  small,  when  R/2as  «  1,  equation  (3)  becomes 

G.^.-rjJjL.  (4) 

For  crystals  with  high  melting  point,  (Tm  -  T0)  is  large  and  the  heat  transfer  coefficient  h 
is  increased  by  radiation  heat  transfer.  The  attainable  Gs  can  be  as  high  as  100°C/cm  or 
more.  For  low  melting  point  crystals,  external  cooling  has  to  be  applied  to  achieve  steep 
gradients. 
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It  should  be  noted  that  the  analytical  solutions  cannot  completely  explain  and  predict  the 
crystal  growth  in  a  practical  vertical  Bridgman  system.  This  is  mainly  56 
because  the  simplifying  assumptions  usually  tend  to  discard  some  important  details 
which  control  crystal  growth.9'10  Accordingly,  numerical  studies  have  been  performed  to 
gain  more  understandings  in  the  vertical  Bridgman  growth  process.7,11'18  Studies  have 
been  focused  on  controlling  the  solid-melt  interface  shape,11,13  models  of  the  heat  and 
mass  transfer, 14,15,18  and  the  influence  of  the  thermal  environment  and  process 
parameters.16'17  Accurate  material  thermal  properties  and  operating  environment  are 
required  for  successful  numerical  predictions. 

1.2  Physical  Vapor  Transportation  Technique 

The  crystal  growth  by  the  PVT  method  is  typically  conducted  in  closed  or  semi-closed 
systems.  If  the  vapor  pressure  of  a  material  at  a  desired  growth  temperature  exceeds 
10'2  torr,  it  may  be  efficiently  grown  by  the  PVT  method.  In  the  PVT  process,  an 
ampoule  containing  a  polycrystalline  source  of  the  desired  material  is  heated  to  a 
temperature  at  which  the  material  starts  to  sublime  at  a  rate  leading  to  crystal  growth. 
The  ampoule  is  typically  placed  in  a  furnace  that  provides  a  temperature  gradient  over 
the  length  of  the  ampoule,  so  that  the  polycrystalline  source  material  sublimes  at  the 
higher  temperature  and  deposits  at  the  end  of  the  ampoule  where  the  temperature  is 
lower.  This  temperature  difference  causes  the  different  vapor  pressure  at  two  ends  of 
the  ampoule  and  forces  the  material  transport  from  the  source  to  the  cooler  end.  The 
PVT  system  can  adopt  horizontal  or  vertical  set  up.  A  typical  horizontal  PVT  system  is 
illustrated  in  Figure  2. 


Figure  2.  Schematic  diagram  of  physical  vapor  transport  system.20 

For  crystal  growth  by  PVT  method,  the  growth  temperature  T2  should  exceed  about 
one-third  of  the  melting  temperature  of  the  crystal.6  The  temperature  Ti  is  fixed  so  that 
the  crystal  grows  at  an  acceptable  rate  which  can  be  deduced  from  simple  theory  as 
follow.  The  vapor  pressure  can  be  expressed  in  the  form  of 
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then  the  rate  of  deposition  in  an  ideal  case  with  an  evacuated  ampoule  is6 


f  =  P0 


-Jf-lp 

2 nRT2  / 


(6) 


where  AHS  is  the  latent  heat  associated  with  the  vaporization  of  a  solid  below  its  melting 
point  (cal/g),  p  is  the  density  of  the  crystal  (g/cm3),  M  is  the  molecular  weight  of  the 
vapor  (g/mole),  and  R  is  the  gas  constant  (cal/(K-mole)).  Equation  (6)  indicates  that  the 
growth  rate  of  PVT  method  can  be  controlled  by  adjusting  the  temperature  difference  at 
two  ends  of  the  ampoule.  Normally  inert  gases  (N2,  Ar)  are  infused  into  the  ampoule 
after  vacuum  treatment  to  assist  the  nucleation  of  PVT  crystal  growth.21"23  Since  PVT 
growth  is  carried  out  at  temperature  much  lower  than  the  melting  point  of  the  material,  it 
has  multiple  benefits  of24:  (a)  reduced  diffusion  of  impurities  from  containers,  (b) 
reduced  defects  such  as  voids  and  dislocation  densities,  and  (c)  effective  separation  of 
impurities  with  low  vapor  pressure  during  transport.  The  major  drawbacks  for  PVT 
methods  are  small  crystal  volume  and  uncontrolled  nucleation.  In  practice,  unless  there 
is  a  large  temperature  gradient  near  the  cooler  end  of  the  ampoule,  multiple  nucleations 
occur  and  the  products  at  the  cooler  end  are  needles,  platelets,  or  polycrystalline 
formed  by  many  small  crystals. 

1.3  Modified  Vertical  Bridgman  Growth 

In  this  study,  Bil3  single  crystals  were  mainly  grown  by  a  modified  vertical  Bridgman 
method.  Pyrex  glass  ampoules  were  selected  as  the  crystal  growth  containers  since  the 
annealing  and  softening  points  of  the  glass  are  560°C  and  815°C  respectively,  which 
are  much  higher  than  the  maximum  temperature  (440°C)  used  in  the  growth  of  Bil3 
single  crystal.  The  Pyrex  glass  tube  was  soaked  in  Decon  Contrad®  70  for  12  hours  to 
effectively  remove  the  organics,  rinsed  10  times  with  high  purity  Dl  water  (p  =  13.6  MQ- 
cm)  and  baked  at  120°C  for  12  hours.  One  end  of  the  cleaned  glass  tube  was  melted 
and  shaped  to  a  sharp  tip  using  a  propane  hand  torch.  Bil3  powder  (99.999%,  MV 
Laboratories)  was  weighed  and  loaded  in  this  tipped  Pyrex  glass  ampoule.  The 
ampoule  needs  to  be  evacuated  and  sealed  to  get  rid  of  the  water  and  oxygen  as  Bil3 
will  react  with  them  at  elevated  temperature.  To  do  that,  the  ampoule  was  first 
vacuumed  to  a  pressure  of  approximately  4  x  10'5  mbar.  Heat  was  then  gently  applied 
on  the  ampoule  wall  using  a  propane  hand  torch  to  help  evaporate  the  water.  This  step 
has  to  be  done  very  carefully  to  avoid  heating  up  the  Bil3  powder.  Ar  gas  was  then  back 
filled  in  the  ampoule  to  -700  mbar.  The  purpose  of  infusing  Ar  gas  is  to  make  sure  any 
residual  gas  after  vacuum  is  inert.  The  cycle  of  pumping  down,  heating  ampoule  wall, 
and  backfilling  Ar  gas  was  repeated  three  times,  the  final  evacuation  took  around 
another  2-3  hours  to  take  the  pressure  down  to  10"6  mbar.  The  ampoule  was  then 
sealed  using  a  propane  hand  torch.  To  help  seal  the  ampoule,  a  neck  area  was 


preformed  in  air  before  the  ampoule  was  pumped  down.  As  the  soft  glass  will  collapse 
inwards  due  to  the  inner  vacuum,  the  neck  area  will  help  the  glass  collapse  into  the 
opposing  wall  without  breaking  the  seal.  After  sealing,  the  ampoule  was  vertically  fixed 
in  a  standing  frame  and  placed  in  a  programmable  24  heating  zone  furnace  (EDG-13, 
Mellen  Company).  The  length  of  each  heating  zone  is  one  inch  and  the  temperature  of 
each  heating  zone  can  be  independently  controlled  by  a  computer  program,  which 
allows  a  precisely  controlled  temperature  profile  along  the  vertical  axis  of  the  furnace. 
The  powder  was  first  heated  and  kept  at  440°C  for  24  hours  then  cooled  to  420°C  at  a 
rate  of  10°C/h.  The  temperature  profile  was  programmed  to  move  up  along  the  vertical 
axis  with  a  fixed  rate.  This  simulated  the  downward  movement  of  the  ampoule  in  the 
conventional  vertical  Bridgman  method.  The  set  up  of  this  modified  Bridgman  system 
provides  two  advantages  over  the  conventional  one.  First,  possible  vibration  during 
crystal  growth  is  eliminated  since  the  ampoule  is  kept  static  in  the  new  system.  Second, 
a  hole  is  no  longer  needed  at  the  top  of  the  furnace  to  provide  the  pass  way  for  the 
hanging  wire.  Thus,  less  thermal  disturbance  is  expected.  As  a  result,  the  modified 
Bridgman  system  can  provide  a  more  stable  crystal  growth  environment  leading  to 
better  crystal  qualities. 

As  discussed  previously,  temperature  gradient  and  growth  rate  are  two  critical 
parameters  for  crystal  growth  by  vertical  Bridgman  technique.  The  best  way  to  start  the 
crystal  growth  is  to  assess  the  relationship  between  the  growth  rate  and  temperature 
gradient  at  solid-liquid  interface  established  in  equation  (1).  However,  it  is  difficult  to  do 
that  in  the  case  of  Bil3  due  to  the  unavailability  of  its  thermal  conductivity  values.  This  is 
probably  because  Bil3  will  dissociate  to  Bi  and  l2  at  250°C  -  300°C,25  and  l2  is  very 
corrosive  with  extremely  high  reactivity  to  metals,  which  could  destroy  the  instrument  for 
the  thermal  conductivity  measurement.  Therefore,  in  this  study,  different  growth  rates 
and  temperature  gradients  for  growing  crystals  with  0.75  inch  diameter  were 
investigated  on  an  educated  initial  guess  basis.  Previous  studies  reported  that  Bil3 
single  crystals  have  been  grown  at  the  conditions  of  1  mm/h  growth  rate  with  15°C/cm 
temperature  gradient1  and  1  mm/h  growth  rate  with  5.7°C/cm  temperature  gradient.2  As 
equation  (4)  indicates  that  the  temperature  gradient  depends  on  the  thermal  properties 
of  the  furnace,  Bil3  single  crystal  grown  in  different  furnaces  could  have  very  different 
temperature  gradient.  The  highest  temperature  gradient  allowed  in  the  EDG  24-zone 
furnace  is  around  18°C/cm.  The  temperature  variation  in  each  zone  is  controlled  within 
±0.5°C,  which  corresponds  to  a  gradient  variation  of  ±0.2°C/cm.  The  investigated 
temperature  gradients  were  then  selected  as  5,  10,  and  15°C/cm  to  make  a  distinct 
difference  between  them.  For  each  temperature  gradient,  the  crystal  growth  has  been 
conducted  at  a  speed  of  0.5,  1 ,  and  2  mm/h. 

1 .4  Physical  Vapor  T ransport  Growth 

For  the  PVT  crystal  growth,  the  Pyrex  glass  ampoule  was  cleaned  using  the  same 
procedure  described  for  the  MVB  growth.  A  charge  of  5  g  Bil3  powder  was  loaded  into  a 
cleaned  ampoule  with  diameter  of  0.75  inches.  The  work  charge  was  first  vacuumed  to 
around  4  x  10'5  mbar  and  cleaned  with  ultra  high  purity  grade  Ar  gas  three  times,  then 
finally  sealed  with  ultra  high  purity  grade  Ar  at  700  mbar.  Following  the  crystal  growth 
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conditions  used  by  Cuna  et  al ,,21  the  ampoule  end  with  Bib  powder  was  placed  at 
310°C  and  the  other  end  was  placed  at  room  temperature.  The  growth  time  was  7  days 
for  5  g  of  powder.  After  the  growth,  the  furnace  was  cooled  down  at  a  rate  of  1°C/min. 
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2.  Crsytal  Growth  Results  and  Discussion 

2.1  MVB  Crystal  Growth 

2.1.1  Determination  of  Temperature  Profile 

The  concept  of  vertical  Bridgman  growth  is  quite  simple.  The  powder  is  first  heated  to 
the  melting  temperature  and  soaked  for  a  period  of  time.  Then  the  crystal  growth  is 
conducted  under  certain  temperature  gradient  and  moving  rate.  Finally  the  entire  growth 
system  is  cooled  back  to  room  temperature.  While  the  temperature  gradient  and  growth 
rate  are  vital  for  producing  high  quality  single  crystals,  the  heat  up  and  cool  down  cycles 
are  also  important  for  a  successful  growth  run. 

As  the  furnace  thermocouples  are  located  at  the  furnace  wall  while  the  growth  ampoule 
is  in  the  center,  the  temperature  difference  is  presented  between  the  zone  temperature 
reading  and  the  actual  temperature  of  the  ampoule.  To  determine  this  difference,  the 
central  temperature  was  recorded  by  a  K-type  thermocouple  placed  along  the  central 
axis  of  the  furnace.  It  was  observed  that  the  central  temperature  was  around  20°C 
higher  than  the  temperature  reading  at  the  furnace  wall.  To  ensure  the  complete  melting 
of  Bib,  the  highest  temperature  for  the  soaking  period  was  set  as  440°C.  This  is  based 
on  the  fact  that  ideally  Bib  will  melt  at  408°C  and  the  actual  soaking  temperature  will  be 
around  460°C. 

The  heating  rate  is  not  a  major  concern  for  Bib  as  the  processing  temperature  is  in  the 
medium  range  and  Pyrex  ampoule  has  a  very  low  thermal  expansion  coefficient  (3.2  x 
10'6/°C).  The  only  concern  is  to  make  sure  that  the  furnace  temperature  does  not  go  up 
too  fast  or  too  slow,  as  both  conditions  require  additional  time  to  reach  thermal 
equilibrium.  Experimentally,  the  heating  rate  has  been  tested  at  1 ,  3,  and  5°C/min.  It 
was  found  that  at  3  and  5°C/min  heating  rates  the  furnace  temperature  could  not  follow 
the  set  temperature.  The  temperature  lags  were  of  around  10  and  30°C  for  the  heating 
rate  of  3  and  5°C/min  respectively.  When  the  heating  rate  was  lowered  to  1°C/min,  the 
furnace  temperature  completely  followed  the  set  point  with  ±1°C  variations.  Based  on 
the  test  result,  the  heating  rate  was  set  at  1°C/min  for  all  the  crystal  growth  runs. 

After  the  furnace  reached  the  maximum  temperature,  the  powder  in  the  growth  ampoule 
may  still  be  lower  than  this  point  since  heat  transfer  takes  time.  A  soak  period  has  to  be 
set  up  to  allow  every  component  inside  the  furnace  reach  the  thermal  equilibrium.  A  K- 
type  thermal  couple  was  attached  at  the  outer  wall  of  the  ampoule  to  test  the  equilibrium 
time.  It  was  observed  that  the  temperature  was  stabilized  in  the  first  30  minutes  after 
reaching  the  maximum  temperature.  However,  the  Bib  melt  inside  the  ampoule  may 
need  longer  time  to  reach  thermal  equilibrium.  A  soaking  time  of  12  hours  will  be  safe 
enough  to  thermally  stabilize  the  whole  growth  system. 

After  soaking,  the  furnace  was  cooled  down  to  410°C  at  a  rate  of  1°C/min.  The 
corresponding  central  temperature  was  around  430°C.  This  is  to  make  sure  the  overall 


furnace  temperature  is  still  above  the  melting  point  of  Bil3.  Starting  from  this  point,  the 
crystal  began  to  grow  at  a  specific  growth  rate  and  temperature  gradient  described 
previously.  The  temperature  at  which  the  crystal  growth  ends  is  another  important 
parameter  needed  to  be  determined.  If  the  crystal  growth  ends  at  a  high  temperature, 
lattice  disorder  may  happen  due  to  the  fast  cooling  of  the  atoms  with  high  kinetic 
energies.  On  the  other  hand,  ending  the  crystal  growth  at  a  very  low  temperature  is  time 
wasting.  It  is  noted  that  the  crystal  growth  of  Bil3  by  PVT  method  happens  at  ~250°C,4 
which  implies  that  atoms  at  this  temperature  barely  have  enough  kinetic  energies  to 
move  around.  Therefore,  in  this  study,  the  crystal  growth  was  set  to  end  at  200°C,  which 
corresponding  to  a  furnace  temperature  of  180°C.  The  crystal  after  growth  was  cooled 
to  room  temperature  at  a  rate  of  1°C/min.  A  temperature  profile  was  created  based  on 
the  above  discussion,  as  illustrated  in  Figure  3.  The  total  growth  time  varies  with  the 
temperature  gradient  and  growth  rate.  For  example,  the  growth  time  for  a  crystal  grown 
at  10°C/cm  and  0.5  mm/h  is  around  20  days. 


Figure  3.  Growth  recipe  for  the  Bil3  crystal  growth  at  10°C/cm  and  0.5  mm/h. 

2.1.2  Crystal  Growth  Result 

For  crystal  growth  in  0.75  inch  diameter  ampoules,  it  was  determined  that  only  three 
different  combinations  of  growth  rates  and  temperature  gradients  produced  large  single 
crystals:  10°C/cm  and  0.5  mm/h;  10°C/cm  and  1  mm/h;  and  15°C/cm  and  0.5  mm/h 
(Table  1).  The  typical  polycrystalline  rods  and  cleavage  planes  of  crystals  obtained  at 
different  growth  conditions  are  shown  in  Figure  4. 


Table  3-1 


Temperature  Gradient  (oC/cm) 

5 

10 

15 

Growth 

0.5 

Potycrystalline 

Single  crystal 

Single  crystal 

Rate 

1 

Pol /crystalline 

Single  crystal 

Polycrystalline 

(mm/h) 

2 

Pol /crystalline 

Polycrystalline 

Pol /crystal  line 

It  is  noticed  that  single  crystals  were  obtained  at  the  same  temperature  gradient 
(10°C/cm)  but  different  growth  rates  (0.5  and  1  mm/h)  or  the  same  growth  rate  (0.5 
mm/h)  but  different  temperature  gradients  (10  and  15°C/cm).  It  is  therefore  clear  that 
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equation  (1)  cannot  explain  these  experimental  results.  However,  one  has  to  keep  in  is 
established  on  the  steady  state  condition.  A  real  crystal  growth  always  happens  at  a 
non-steady  state,  although  it  may  be  very  close  to  a  steady  state.  Therefore,  the 
requirement  for  crystal  growth  becomes 


psLR  <  KSGS  - KlGl  (7) 

which  suggests  as  long  as  the  overall  heat  flux  flowing  away  from  the  growth  interface 
exceeds  the  rate  of  latent  heat  evolution,  the  growth  of  single  crystal  could  happen.6 
This  explains  why  Bil3  single  crystals  were  obtained  at  the  three  conditions  listed  in 
Table  1. 


Figure  4.  Single  crystals  produced  by  MVB  method,  (a)  as  grown  polycrystalline  rod;  (b)  crystal  grown  at 
10°C/cm  and  0.5  mm/h;  (c)  crystal  grown  at  10°C/cm  and  1.0  mm/h;  crystal  grown  at  15°C/cm  and  0.5 
mm/h;  (e)  crystal  grown  at  15°C/cm  and  1.0  mm/h;  (f)  polycrystalline  cross  section  obtained  at  rest 

conditions. 

Based  on  the  above  discussion,  single  crystals  should  also  be  grown  at  15°C/cm  and  1 
mm/h  given  that  single  crystal  was  obtained  at  10°C/cm  and  1  mm/h.  Actually,  the 
solidification  at  this  condition  was  characterized  by  parallel  growth  of  multiple  crystals, 
which  had  a  non-smooth  growth  interface,  as  shown  in  Figure  3-7e.  A  possible  reason 
for  this  kind  of  interface  morphology  was  the  thermal  fluctuations  at  the  growth  interface, 
a  phenomenon  that  has  been  commonly  observed  for  other  single  crystal  growths.26'28 
The  growth  furnace  was  observed  to  have  a  temperature  fluctuation  of  around  ±1°C  at  a 
temperature  gradient  of  15°C/cm,  compared  to  ±0.5°C  at  a  temperature  gradient  of 
10°C/cm.  A  higher  growth  rate  would  cause  more  fluctuations  and  non-smooth  interface 
compared  to  a  lower  growth  rate.  Therefore,  crystal  grown  at  15°C/cm  and  1  mm/h  was 
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under  a  condition  with  high  thermal  fluctuations,  which  caused  the  formation  of  a  non¬ 
smooth  interface.  By  decreasing  the  temperature  gradient  to  10°C/cm  or  the  growth  rate 
to  0.5  mm/h,  the  thermal  fluctuations  were  suppressed  and  smooth  growth  face  was 
obtained.  Based  on  the  experimental  results,  the  critical  conditions  to  grow  Bil3  single 
crystals  with  a  diameter  of  0.75  inch  can  be  concluded  as  a  temperature  gradient 
between  10°C/cm  and  15°C/cm,  and  a  growth  rate  between  0.5  mm/h  and  1  mm/h.  A 
single  crystal  with  a  dimension  of  18  x  13  x  5  mm3  has  been  successfully  grown  at  the 
temperature  gradient  of  10°C/cm  and  growth  rate  of  0.5  mm/h,  as  shown  in  Figure  5. 


Figure  5.  Bil3  single  crystal  with  dimension  Of  18  x  13  x  5  mm3,  the  largest  crystal  size  ever  reported. 

This  is  the  largest  crystal  size  ever  reported  for  Bil3.  The  previously  reported  largest  Bil3 
single  crystal  had  a  dimension  of  12  x  12  x  4  mm3.29 

2.1.3  PVT  Crystal  Growth 

The  growth  ampoule  and  a  single  crystal  platelet  obtained  after  PVT  growth  are  shown 
in  Figure  6.  Many  tiny  crystals,  needles,  and  platelets  were  found  to  stack  up  together  in 
the  ampoule  where  the  temperature  was  measured  to  be  around  230°C,  which  was 
determined  to  be  the  growth  temperature  of  the  single  crystal  platelets.  The  platelets 
were  carefully  separated  from  the  polycrystalline  mixture  by  using  tweezers  and  a  razor 
blade.  The  typical  platelet  has  a  surface  area  of  5  -  20  mm2  and  a  thickness  of  100  - 

200  pm. 


Figure  6.  Bil3  single  crystals  grown  by  the  PVT  method. 
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3.0  Crystal  Characterization 

A  polycrystalline  ingot  was  obtained  after  the  crystal  growth  by  the  MVB  method.  Single 
crystal  samples  were  cut  from  the  ingot  along  the  cleavage  plane  (001)  with  various 
cross  sectional  areas  and  thicknesses.  For  crystal  grown  by  the  PVT  method,  single 
crystal  platelets  with  various  dimensions  were  carefully  collected.  Single  crystal  samples 
obtained  from  both  methods  were  characterized  by  the  following  techniques. 

3.1  Characterization  Techniques 

3.1.1  X-ray  Diffraction 

X-ray  diffraction  using  CuKa  radiation  (Philips  APD  3720)  was  used  to  determine  the 
crystallinity  and  phase  purity  of  Bil3  single  crystals.  The  power  of  the  X-ray  generator 
was  set  to  40  kV  and  20  mA.  The  20  angle  was  scanned  from  10°  to  60°  with  a  step  size 
of  0.02° 

3.1.2  Vickers  Hardness 

The  Vickers  hardness  of  the  grown  single  crystals  was  measured  by  a  micro  hardness 
tester  (Buehler  MicroMet  II).  Each  sample  was  measured  for  5  indentations  and  the 
applying  force  was  50  g. 

3.1.3  Etch  Pit  Density 

The  dislocation  density  of  the  crystal  surface  was  examined  by  EPD  measurement 
using  a  Kl  solution  as  the  etchant.  To  determine  the  ideal  etchant  concentration  for  EPD 
measurement,  single  crystal  samples  were  etched  in  5%,  10%,  15%,  and  20%  weight 
percent  of  Kl  solution  respectively,  and  the  etching  time  varied  from  10  sec  to  150  sec. 
After  etching,  the  sample  was  rinsed  with  acetone  and  dried  at  120°C  for  10  min.  The 
etched  sample  was  weighed  to  determine  the  weight  loss  per  unit  surface  area.  The 
time  dependence  of  the  weight  loss  was  calculated  for  each  etchant  concentration.  The 
ideal  etchant  concentration  was  determined  at  which  the  etchant  can  effectively  remove 
the  mass  of  Bil3  single  crystal.  After  that,  Bil3  single  crystals  grown  by  the  MVB  and  PVT 
methods  were  etched  for  a  specific  time.  The  density  of  pits  was  examined  by  an  optical 
microscope  (Nikon  Eclipse  LV100).  For  each  sample,  five  randomly  selected  areas 
were  studied  to  count  the  EPD. 

3.1.4  X-ray  Rocking  Curve 

The  quality  of  the  grown  crystals  was  inspected  by  X-ray  rocking  curve  on  a  Philips 
analytical  MRD  X'Pert  System,  as  shown  in  Figure  7.  The  presence  of  defects  such  as 
dislocations  and  small  angle  grain  boundaries  will  cause  local  variation  of  lattice  spacing 
or  plane  orientation.  For  X-ray  diffraction  in  the  vicinity  of  a  defect,  Bragg’s  law  cannot 
be  simultaneously  satisfied  by  the  perfect  lattice  and  the  distorted  one,  which  causes 
the  broadening  of  X-ray  diffraction  peak.  By  performing  X-ray  rocking  curve  analysis, 
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the  peak  broadening  of  a  specific  lattice  plane  can  be  measured  by  scanning  the 
incident  angle  w  at  a  fixed  receiving  angle  20,  as  illustrated  in  Figure  8.  The  full  width  at 
half  maximum  (FWHM)  of  this  peak  is  generally  used  to  evaluate  the  quality  of  a  single 
crystal.30'32  The  rocking  curve  ou-scan  was  conducted  at  a  rate  of  0.002°/s  with  a  step 
size  of  0.001°.  Each  sample  was  scanned  twice  by  rotating  the  sample  90°  (<i>  angle  in 
Figure  8)  relative  to  the  previous  scan,  as  illustrated  in  Figure  9. 


Figure  7.  X’Pert  system  for  X-ray  rocking  curve  measurements. 
Diffracted  X-ray  \  <j>  Incident  X-ray 


Figure  8.  Illustration  of  rocking  curve  w-scan. 


Figure  9.  Illustration  of  the  sampling  for  rocking  curve  w-scan.  The  shadowed  bands  numbered  1  and  2 

represent  the  x-ray  incident  area  for  each  scan. 
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3.1.5  Scanning  Electron  Microscopy  (SEM) 

The  surface  microstructure  of  the  as  grown  Bil3  single  crystals  was  characterized  using 
SEM  (JEOL  6335F).  The  freshly  cleaved  surface  was  sputtered  with  carbon  film  with 
thickness  of  around  20  nm.  The  SEM  was  operated  with  an  accelerating  voltage  of  15 
kV,  a  probe  current  of  10  pA,  and  a  working  distance  of  approximately  15  mm. 

3.1.6  Atomic  Force  Microscopy  (AFM) 

To  study  the  surface  morphology  of  the  metal  contacts,  the  sample  surfaces  were 
examined  by  AFM  (Asylum  Research  AFM).  The  AFM  tip  from  Veeco  has  a  triangular 
cantilever  and  a  nominal  tip  radius  of  20  nm.  The  normal  stiffness  of  the  tip  cantilever  is 
0.58  N/m.  The  topography  images  were  acquired  in  contact  mode  at  a  fixed  applied 
load  of  several  nano  Newton. 

3.1.7  X-ray  Photoelectron  Spectroscopy  (XPS) 

The  stoichiometry  of  the  grown  Bil3  single  crystal  is  determined  by  XPS  spectrometer 
using  an  Omicron  XPS  system.  The  XPS  spectra  were  collected  from  a  freshly  cleaved 
Bil3  crystal  surface  with  an  Al  Ka  monochromatic  X-ray  (1486.7  eV).  The  spectra  were 
recorded  at  a  sampling  area  of  1 .3  mm2  and  a  take-off  angle  of  55°  with  pass  energy  of 
22  eV  (step  size  0.05  eV,  integration  time  1  sec).  The  chamber  pressure  was  kept 
below  5  x  10'10  mbar  during  measurement.  The  instrument  resolution  was  around  0.1 
eV.  The  peak  energies  were  calibrated  using  a  Ag  3d5/2  peak  at  368  eV.  The 
composition  of  bismuth  and  iodine  was  estimated  by  comparing  the  photoelectron  peak 
intensities  of  the  Bi  4f  with  the  I  3d  core  level  peaks. 

3.1.8  X-ray  Fluorescence  Spectroscopy  (XRF) 

The  stoichiometry  of  the  grown  single  crystal  was  also  determined  by  XRF  using  an 
Orbis  Micro-XRF  spectrometer  (EDAX  Inc.).  The  XRF  spectrometer  was  operated  using 
a  Rh  tube  at  30  kV.  The  spot  size  of  the  X-ray  was  30  pm.  The  composition  of  bismuth 
and  iodine  was  estimated  by  comparing  the  X-ray  intensities  of  the  Bi  M-line  with  the  I 
L-line  in  the  measured  spectrum. 

3.2  Results  and  Discussion 

3.2.1  Crystal  Structure  and  Phase  Purity 

X-ray  diffraction  was  used  to  determine  the  crystal  structure  and  phase  purity  of  single 
crystals  grown  by  the  two  techniques  described  before.  All  the  diffraction  patterns  were 
almost  identical  and  a  typical  XRD  pattern  is  shown  in  Figure  10.  Only  003,  006,  009, 
and  0012  reflections  parallel  to  the  (001)  plane  were  observed,  and  the  peak  positions 
matched  well  with  the  published  Bil3  powder  diffraction  data  for  a  rhombohedral  crystal 
structure  (JCPDS  PDF#  48-1 795).29 
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The  lattice  parameters  of  the  single  crystal  grown  by  the  MVB  method  were  determined 
by  performing  20  scans  on  several  selected  planes.  The  scan  was  conducted  by 
rotating  and  tilting  the  single  crystal  to  satisfy  the  Bragg’s  condition  for  a  specific  plane 
that  is  not  parallel  to  the  crystal  surface.  The  planes  and  their  corresponding  20  angle 
are  listed  in  Table  3-2.  The  lattice  parameters  were  calculated  by  least  squares  fit  with 
values  of  a  =  7.51 6±0. 002  A  and  c  =  20.697±0.003  A,  which  matched  well  with  the 
published  powder  diffraction  data  of  a  =  7.51 9  A  and  c  =  20.721  A.1 0  The  XRD  results 
verified  that  Bil3  single  crystals  have  been  successfully  grown  by  the  MVB  and  PVT 
methods. 


Figure  10.  A  typical  XRD  pattern  for  Bil3  single  crystals. 


Table  2.  20  Angles  of  Selected  Planes. 


Planes  (hkl) 

003 

104 

107 

116 

018 

20 

12.790 

21.920 

33.258 

35.205 

37.371 

3.2.2  Vickers  Hardness 

Vickers  hardness  test  is  a  method  to  determine  the  hardness  of  a  material  by  measuring 
the  size  of  an  indentation  produced  under  force  by  a  pyramid  shape  diamond  indenter.33 
The  angle  between  the  two  opposite  faces  of  the  pyramid  is  136°.  After  removal  of  the 
load,  the  two  diagonals  (dl,  d2)  of  the  indentation  on  the  surface  of  the  material  are 
measured  under  a  microscope  and  their  average  is  calculated  (Figure  11).  The  area  of 
the  sloping  surface  of  the  indentation  is  calculated  by 

,  d 2  d2  ^ 

A  = - f - t^i  « -  (°) 

2sin(l36^)  1-8544 


where  d  is  the  mean  diagonal  of  the  indentation  with  a  unit  of  mm.  The  Vickers 
hardness  (HV  number)  is  then  determined  by 
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HV  =  — 
A 


F  ~  1.8544F 
A  ~  d2 


(9) 


where  F  is  the  applied  force  with  a  unit  of  Kg. 


Force 


l 


Figure  11.  Illustration  of  the  Vicker’s  hardness  test. 


The  Vickers  hardness  on  the  (001)  surface  of  single  crystals  grown  by  the  MVB  method 
was  measured  and  is  shown  in  Table  3.  Single  crystals  grown  by  the  PVT  method  were 
not  subjected  to  the  hardness  test  as  the  crystals  are  too  thin.  The  average  HV  numbers 
of  the  single  crystals  grown  at  different  conditions  showed  no  significant  difference.  The 
results  indicated  that  Bil3  single  crystal  is  soft  with  a  HV  number  between  12  and  15.  As 
a  comparison,  the  HV  number  of  a  Hgl2  (001)  surface  is  ~22,34  a  316L  stainless  steel  is 
140. 35 


Table  3.  Vickers  hardness  of  (001)  surface  of  Bil3  single  crystals  grown  by  the  MVB  method. 


Meas. 

10°C/cm,  0.5  mm/h 

dl  d2  HV# 

10°C/cm,  1  mm/h 

dl  d2  HV# 

15°C/cm.  0.5  mm/h 
dl  d2  HV# 

1 

0.088 

0.089 

12 

0.082 

0.078 

14 

0.077 

0.079 

15 

2 

0.073 

0.088 

14 

0.093 

0.085 

12 

0.083 

0.080 

14 

3 

0.082 

0.089 

13 

0.088 

0.085 

12 

0.086 

0.079 

14 

4 

0.091 

0.088 

12 

0.082 

0.080 

14 

0.075 

0.080 

15 

5 

0.086 

0.087 

12 

0.082 

0.087 

13 

0.084 

0.080 

14 

Aver. 

13±0.15 

13±0.16 

14±0.18 

3.2.3  Etch  Pit  Density  Measurement 

The  EPD  measurement  is  generally  utilized  to  determine  the  quality  of  semiconductor 
single  crystals.36  38  An  etching  solution  is  applied  on  the  surface  of  a  crystal  and  the 
etching  rate  increases  at  dislocation  sites  resulting  in  the  formation  of  pits.  Therefore, 
EPD  is  often  used  to  evaluate  the  dislocation  density  of  single  crystals.39  In  this  study, 
the  etchant  was  Kl  solution  as  it  was  routinely  used  to  reveal  dislocations  in  Hgl222  and 
Pbl240.  The  time  dependence  of  the  weight  loss  for  different  etchant  concentrations  is 
shown  in  Figure  12. 
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Figure  12.  Bil3  etched  in  different  concentrations  of  Kl  solution.  The  unit  of  weight  loss  is  mg/cm2. 

It  was  observed  that  the  initial  etching  rate  increased  with  the  etchant  concentration. 
The  weight  loss  reached  a  maximum  after  etching  approximately  2  minutes  for  all 
solutions.  The  maximum  weight  loss  was  around  27  mg/cm2,  which  corresponded  to  an 
etched  depth  of  approximately  50  pm  for  Bib.  The  weight  loss  curve  showed  almost  no 
difference  once  the  etchant  concentration  was  higher  than  15%.  Therefore,  15%  Kl 
solution  was  selected  as  the  etchant  for  the  EPD  study. 


The  optical  images  of  Bil3  single  crystal  surfaces  after  etched  in  15%  Kl  solution  are 
shown  in  Figure  13.  The  pits  had  a  curved  triangular  shape,  which  represented  the 
atomic  arrangement  of  (001)  plane.  The  size  of  the  pits  was  around  60  pm  (the  side 
length  of  the  triangle)  after  etched  for  1 0  sec.  It  was  observed  that  the  size  of  the  pits 
increased  with  the  increasing  etching  time.  After  etched  for  30  sec,  the  pits  started  to 
overlap  with  each  other  and  the  number  of  pits  became  uncountable.  Therefore  the 
etching  time  for  EPD  measurements  was  determined  to  be  10  and  20  sec. 
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Figure  1 3.  Optical  images  of  the  Bil3  single  crystal  surface  after  etched  for  time  of  (a)  1 0  sec;  (b)  20  sec; 

(c)  30  sec;  (d)  40  sec;  (e)  50  sec;  (f)  60  sec. 

The  typical  images  of  etch  pits  for  BN3  single  crystals  are  shown  in  Figure  14.  No 
obvious  etched  pits  were  observed  for  crystals  grown  by  the  PVT  method  even  when 
the  etching  time  was  extended  from  10  to  60  sec.  This  indicated  that  crystals  grown  by 
the  PVT  method  had  much  lower  dislocation  density  than  crystals  grown  by  the  MVB 
method. 
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Figure  14.  Etch  pits  of  Bil3  single  crystals  grown  by  MVB  and  PVT  methods,  (a)-(c)  MVB  crystals  grown  at 
different  conditions  and  etched  for  10  sec;  (d)  PVT  crystal  etched  for  10  sec;  (e)  PVT  crystal  etched  for  60 

sec. 

When  counting  the  number  of  the  pits,  if  the  pit  crosses  the  corner  of  the  area,  it  is 
counted  as  1/4  pit;  if  the  pit  crosses  the  side,  it  is  counted  as  1/2  pit.  The  pit  number 
was  counted  in  an  area  of  6.3  x  10-3  cm2.  The  results  of  EPD  measurements  are 
summarized  in  Table  4. 

Table  4.  Etch  pit  density  of  single  crystals  grown  by  MVB  method. 


Growth 

condition 

Etching 
lime  (s) 

Area  1 

Number  of  pits 

Area  2  Area  3  Area  4 

Area  5 

Average 

pit# 

EPD 

(104  cm-2) 

10°C/cm 

10 

65.50 

60.25 

81.50 

77.00 

68.50 

70.55 

1 .1  ±0.1 

Q.5nm/h 

20 

83.00 

72.50 

76.00 

76.50 

82.25 

78.45 

1.2±0.1 

10°C/cm 

10 

87.25 

73.50 

105.50 

88.00 

95.50 

89.95 

1 .4  ±0.2 

Imm/h 

20 

91.50 

62.50 

87.00 

97.25 

85.50 

88.75 

1.4±0.1 

1 5°C/cm 

10 

85.00 

90.25 

05.50 

71.25 

83.50 

83.1 

1.3±0.1 

D.5nm/h 

20 

102.50 

86.00 

93.50 

88.00 

77.50 

89.5 

1.4±0.1 

It  was  observed  that  the  number  of  pits  varied  from  one  area  to  another,  which  indicated 
the  inhomogeneous  distribution  of  dislocations  in  the  MVB  crystal.  The  EPD  didn’t  show 
an  apparent  increase  for  a  longer  etching  time  (from  10  sec  to  20  sec)  and  the  extended 
time  mainly  increased  the  size  of  the  pits,  which  is  more  obviously  shown  in  Figure  13. 
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The  EPD  of  Bil3  single  crystals  grown  by  the  MVB  method  was  on  the  order  of  104  cm-2 
and  the  crystal  grown  at  10°C/cm  and  0.5  mm/h  showed  the  smallest  EPD.  The  EPD 
value  of  Bil3  is  comparable  with  other  semiconductor  crystals  such  as  Hgl2  (~103  - 105 
cm'2)41'43  and  CZT  (~104  -  105  cm'2).4445 

3.2.4  X-ray  Rocking  Curve 

X-ray  rocking  curve  characterization  is  a  nondestructive  analysis  method  compared  to 
the  EPD  measurement.  In  this  study,  the  X-ray  rocking  curves  of  Bil3  (006)  plane  were 
measured  as  this  plane  showed  the  highest  intensity  in  X-ray  diffraction  (Figure  10).  The 
X-ray  rocking  curve  w-scans  of  Bil3  single  crystals  grown  by  MVB  and  PVT  methods  are 
shown  in  Figure  15.  Features  of  extra  peaks  and  shoulders  were  observed  in  all  the 
crystals,  which  suggested  that  crystal  planes  in  all  these  samples  were  not  perfectly 
aligned.  Defects  such  as  dislocations  and  twin  planes  may  account  for  these  features. 
Multiple  small  crystals  with  orientation  close  to  [001]  may  also  exist  in  these  crystals. 

The  plot  with  log  scale  can  reveal  the  low  intensity  features  more  clearly.  The  rocking 
curves  scanned  at  two  perpendicular  directions  for  the  same  sample  showed  different 
shape  and  peak  width.  This  indicated  that  the  distribution  of  the  defects  was 
inhomogeneous  in  these  crystals. 
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Figure  15.  X-ray  rocking  curve  w-scans  of  the  single  crystals,  (a)  -  (c):  single  crystals  grown  by  MVB  at 
(a)  10°C/cm  and  1  mm/h;  (b)  10°C/cm  and  0.5  mm/h;  and  (c)  15°C/cm  and  0.5  mm/h.  (d)  Single  crystal 
grown  by  PVT  method.  Solid  line  -  linear  intensity;  circles  -  log  scale  intensity. 


The  FWHM  of  the  rocking  curves  for  each  sample  is  summarized  in  Table  5.  The 
crystals  grown  by  MVB  method  showed  similar  rocking  curve  FWHMs  thus  the  crystal 
quality  could  not  be  differentiated  by  this  technique.  On  the  other  hand,  the  crystal 
grown  by  PVT  method  showed  overall  smoother  peaks  and  smaller  peak  width 
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compared  to  the  crystals  grown  by  MVB  method.  This  suggested  that  crystals  grown  by 
PVT  method  have  better  quality,  which  was  in  agreement  with  the  EPD  result. 


Table  5  FWHM  of  X-ray  rocking  curves  for  Bil3  single  crystals 


MVB  Crystals 

PVT  Crystal 

10°C^crr  and 

1  CPOcm  and  1 

153C/cm  and 

0.5  mn/h 

nm/h 

0.5  mmi'h 

0°  scan 

0.1618 

03438 

0.2292 

0.0415 

90°  scan 

0.1065 

01670 

0.1839 

0.1181 

The  rocking  curve  to-20  scan  was  performed  on  the  MVB  and  PVT  crystals  to  check  the 
distribution  of  the  (006)  d-spacing,  as  shown  in  Figure  16.  It  was  observed  that  the  d- 
spacing  distribution  for  both  crystals  had  wide  ranges,  which  could  partially  account  for 
the  broadening  of  the  rocking  curves. 


Figure  16.  Rocking  curve  iu-20  scan  of  Bil3  single  crystals  indicated  the  wide  distribution  of  d-spacing. 

The  poor  crystalline  perfection  observed  in  Bil3  single  crystals  was  also  commonly 
observed  in  other  iodine  compound  crystals  such  as  Hg I2  and  Pbh  46  This  was  possibly 
due  to  the  softness  of  all  these  materials.  Strains  induced  from  the  thermal  stress  during 
the  cooling  of  the  crystal  could  easily  cause  the  lattice  distortion  thus  the  broadening  of 
the  rocking  curves.  For  Bil3  single  crystal,  similar  FWHM  values  were  also  reported  by 
Nason  and  Keller  who  grew  the  crystals  on  a  seed  using  the  PVT  method.  Their  rocking 
curve  measurement  on  (0012)  plane  showed  FWHM  values  of  600  and  900  sec  of  arc, 
corresponding  to  0.17°  and  0.25°  respectively.29 

The  rocking  curves  of  Bil3  have  exceedingly  large  FWHM  values  compared  to  other 
single  crystals  with  similar  EPD  values.  For  example,  a  typical  CZT  crystal  with  -104  cm' 
2  EPD  value  exhibited  a  rocking  curve  FWHM  of  -0.01° 44  A  high  quality  GaAs  crystal 
with  EPD  less  than  103  cm'2  showed  a  rocking  curve  FWHM  of -0.001  °.47  It  should  be 
noted  that  due  to  the  softness  of  Bil3,  the  crystal  could  also  be  deformed  during  cutting 
and  handling,  as  evidenced  in  Figure  18.  It  is  possible  that  the  peak  broadening  was 
exaggerated  by  these  additional  strains. 
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Figure  18.  Deformation  of  Bil3  caused  by  cutting. 

3.2.5  SEM  Characterization 


The  (001)  surface  morphologies  of  Bil3  single  crystal  were  observed  by  SEM.  The 
surfaces  were  clean,  smooth,  and  free  of  cracks,  as  shown  in  Figure  19.  The 
deformation  caused  by  cutting  was  clearly  observed  at  the  edge  of  the  single  crystal 
(Figure  19a).  The  layered  structure  of  (001)  planes  was  revealed  in  Figure  19b.  The 
grey  square  box  observed  in  Figure  19b  was  the  evidence  of  surface  damage  caused 
by  the  strong  electron  beams,  which  illustrated  the  vulnerability  of  Bil3  to  electrons. 


Figure  19.  SEM  images  of  Bil3  (001)  surfaces  show  (a)  deformation  induced  by  cutting;  (b)  layered 

structure  of  (001)  planes. 


3.2.6  AFM  Characterization 

The  topographic  morphologies  of  the  cleaved  Bil3  surface  were  characterized  by  AFM. 
The  surface  was  observed  to  be  atomically  rough  with  an  average  surface  roughness  of 
0.484  nm  (Figure  20). 
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Figure  20.  AFM  images  of  Bil3  (001)  cleaved  surfaces. 

The  atomic  layers  were  revealed  in  the  topographic  images  with  irregular  edges,  which 
was  possibly  due  to  the  weak  van  der  waal’s  forces  between  the  layers  along  [001] 
direction. 

3.2.7  XPS  Characterization 

The  elemental  composition  near  the  surface  of  Bib  single  crystal  was  characterized  by 
XPS.  The  full  scanned  XPS  spectra  are  shown  in  Figure  21.  On  the  freshly  cleaved 
surface  (1  day  exposure  in  air),  only  carbon  contamination  was  observed.  For  the 
surface  with  two  weeks  exposure  in  atmosphere,  O  Is  peak  was  clearly  identified.  The 
result  suggested  that  Bil3  was  stable  in  air  for  a  short  time,  possibly  several  days.  It  will 
be  eventually  oxidized  after  a  long  time  exposure  in  air.  Therefore,  Bil3  single  crystals 
were  stored  in  inert  gas  environment  to  avoid  the  oxidization. 


Figure  21 .  XPS  spectrum  of  a  freshly  cleaved  Bil3  single  crystal. 

The  Bi  4f  and  I  3d  spectra  of  the  freshly  cleaved  surface  were  collected  with  higher 
resolution,  as  shown  in  Figure  22.  The  The  I/Bi  atomic  ratio  was  calculated  by  the 
intensity  ratio  between  the  I  3d5/2  peak  the  Bi  4f  peaks.  The  calculation  indicated  an 
iodine  deficiency  with  I/Bi  ratio  of  2.78. 
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Figure  22.  XPS  spectra  of  Bi  4f  and  I  3d  for  a  freshly  cleaved  Bil3  crystal  surface. 

3.2.8  XRF  Characterization 

The  XRF  spectrum  of  Bib  is  shown  in  Figure  23.  Quantitative  analysis  was  performed 
by  comparing  the  intensities  of  Bi  M-line  to  I  L-line.  The  result  indicated  an  I/Bi  ratio  of 
2.54.  Considering  XRF  is  not  a  very  sensitive  analytical  method,  an  uncertainty  of -10% 
is  generally  expected  for  the  quantitative  analysis.48'49  The  stoichiometry  determined  by 
XRF  was  consistent  with  the  XPS  characterization.  Since  XRF  has  a  much  higher 
sampling  depth  than  XPS,  the  iodine  deficiency  is  considered  as  the  bulk  defect  of  the 
single  crystal. 


Energy  (keV) 


Figure  23.  XRF  spectrum  of  Bil3  on  a  freshly  cleaved  crystal  surface. 

The  stoichiometric  deviation  of  Bil3  single  crystal  was  also  observed  by  Matsumoto  et  al. 
who  reported  an  I/Bi  ratio  of -2.9  using  XRF  analysis.2  The  iodine  deficiency  issue  has 
been  commonly  observed  in  other  iodide  semiconductor  compound  such  as  Hgl250  and 
Pbl2.51  It  was  also  reported  that  Bil3  tended  toward  dissociation  into  Bi  and  l2  at  250  - 
300°C,5  which  implied  that  this  stoichiometric  deviation  may  come  from  the  evaporation 
of  iodine  after  crystal  growth. 
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3.2.9  Determination  of  Band  Structure 

Bil3  has  been  reported  with  different  band  gap  characteristics  and  values.  The 
experimentally  and  computationally  determined  band  gap  value  of  Bil3  was  ranging  from 
1 .59  eV  to  2.2  eV  as  summarized  in  Table  6.  For  instance,  Vashchenko  and  Timofeev 
measured  the  absorption  spectra  of  Bil3  at  room  temperature  and  reported  an  indirect 
band  gap  of  1 .730  ±  0.005  eV.52'53  Jellison  et  al.  used  two-modulator  generalized 
ellipsometry  to  determine  ordinary  and  extraordinary  band  gaps  of  1 .991  ±  0.005  eV 
and  1 .997  ±  0.002  eV  at  room  temperature.54  Density  functional  theory  (DFT) 
calculations  by  Schluter  et  al.  yielded  a  direct  band  gap  of  ~2.2  eV.55  However,  more 
DFT  calculations  by  Yorikawa  and  Muramatsu  suggested  that  Bil3  has  an  indirect  band 
gap  of  approximately  1 .59  eV  with  a  rhombohedral  crystal  structure.56  The  much  smaller 
band  gap  value  was  ascribed  to  the  well-known,  endemic  underestimation  of  band  gaps 
by  local  density  approximation  (LDA)  calculations.  There  thus  appears  to  be  a 
discrepancy  in  the  electronic  structure  and  band  gap  values  reported  for  Bil3  in  the 
literature. 


Table  6.  Reported  Bil3  Bandgap  Values. 


Crystal 

growth 

method 

Character  izati 
on  method 

Band  gap 
naiture 

Band 
gap  (eV) 

T  (K) 

Theory 

EPM 

Direct 

22 

0 

Theory 

NCPP-LDA 

Indirect 

1.59 

0 

1.730 

293 

Vertical 

Optical 

Indirect 

1.922 

85 

Bridgman 

absorption 

2.029 

20.4 

Vertical 

Optical 

Direct 

1.72 

298 

Bridgman 

absorption 

Vapor 

Optical 

1  73 

298 

transportation 

transmission 

Vapor 

transportation 

Ellipsometry 

1.991 

298 

1.97* 

286 

Vapor 

Optical 

Indirect 

2.072* 

77 

transportation 

reflection 

2.08* 

6 

'value  obtained  as  peak  energy  at  fundamental  absorption  edge 


In  this  study,  the  band  gap  at  room  temperature  was  characterized  by  measuring  optical 
transmission  and  reflection  spectra  at  the  wavelength  of  300  to  800  nm.  The  sample  is  a 
freshly  cleaved  single  crystal  with  thickness  of  0.9  mm.  The  light  source  was  a  Newport 
ARC  lamp  housing  66902.  A  parallel  single  wavelength  light  beam  was  generated 
through  a  CornerstoneTM  260,  Newport  monochromator  and  incident  on  the  crystal 
surface  at  an  angle  of  3°  to  the  crystal  optic  axis.  The  intensity  of  the  incident,  reflected 
and  transmitted  light  was  measured  separately  by  a  UV-silicon  detector  (Newport). 
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3. 2. 9.1  Band  Gap  Determined  by  Optical  Spectra  at  Room  Temperature 

The  UV-Vis  transmission  and  reflection  spectra  were  measured  from  300  nm  to  800  nm 
at  room  temperature,  and  are  shown  in  Figure  24.  The  absorption  coefficient  a  was 
determined  from  the  experimental  spectra  by  solving  the  following  equation:56 


r_(l  -Rfe-ad 
1  -R2e~2ad 


(10) 


where  T  is  transmittance,  R  is  reflectance,  and  d  is  the  sample  thickness.  The  relation 
between  the  absorption  coefficient  a  and  incident  photon  energy  (hv)  can  be  written 

as:57 


(ahv)^  =  A(nv-  Eg)  (11) 

where  a  is  the  optical  absorption  coefficient,  h  is  Planck’s  constant,  v  is  frequency  of  the 
incident  photon,  and  Eg  is  the  optical  band  gap.  The  exponent  n  depends  on  the  type  of 
optical  transition.  For  direct  transition,  n  =  1/2;  for  indirect  transition,  n  =  2.  In  order  to 
determine  the  nature  of  the  band  gap,  both  direct  and  indirect  band  gap  relations  were 
examined. 
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Figure  24.  UV-Vis  reflection  and  transmission  spectra  of  Bil3. 

By  plotting  (ahv)2  or  (ahv)1/2  versus  hv  and  extrapolating  the  linear  part  to  (ahv)2  =0  or 
(ahv)1/2  =  0,  the  optical  band  gap  was  determined,  as  shown  in  Figure  25. 
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Figure  25.  Behavior  of  (ahv)2  or  (ahv)1/2  versus  photon  energy. 

A  stepped  structure  was  observed  in  the  (ahv)1/2  versus  hv  plot,  indicating  an  indirect 
band  characteristic  of  Bil3.  The  band  gap  value  was  determined  by  the  average  of  the 
two  interceptions  at  (ahv)1/2  =  0,  which  gave  an  indirect  band  gap  of  1 .68  ±  0.09  eV  at 
room  temperature.  This  value  is  in  agreement  with  experimental  value  reported  in 
references  50,58'60,  but  smaller  than  the  values  reported  in  references  51,61 ,  as  listed  in 
Table  6. 

3. 2. 9. 2  Effect  of  Temperature  on  Band  Gap 

The  effect  of  temperature  on  the  band  gap  was  determined  by  measuring  UV-Vis 
transmission  spectra  from  298  to  330  K;  the  results  are  shown  in  Figure  26.  It  was 
observed  that  the  absorption  edge  shifts  to  the  lower  energy  side  when  the  temperature 
is  increased.  The  absorption  coefficient  a  was  calculated  assuming  that  the  energy  loss 
of  the  transmitted  light  was  solely  due  to  the  sample  absorption, 

a  =  — -ln(r)  (12) 

d 

where  d  is  the  sample  thickness  and  T  is  transmittance.  The  dependence  of  (ahv)1/2  and 
(ahv)2  on  the  photon  energy  is  plotted  in  Figure  27. 


27 


Wavelength  (nni) 


Figure  26.  UV-Vis  transmission  spectra  of  a  Bil3  single  crystal  at  different  temperatures. 
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Figure  27.  (a)  Dependence  of  indirect  band  gap  model  on  the  photon  energy;  (b)  dependence  of  the 

direct  bandgap  model  on  the  photon  energy. 
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Continuous  curves  instead  of  stepped  characteristics  are  observed  in  the  (ahv)1/2  versus 
hv  plot,  which  indicate  that  the  gap  is  not  indirect  band.  However,  it  is  important  to  note 
that  compared  with  the  optical  spectra,  here  the  reflectance  was  not  taken  into  account 
for  the  calculation  of  absorption  coefficient.  This  is  particularly  relevant  because 
equation  1 1  is  only  valid  for  a  sample  with  negligible  reflectivity,  which  is  not  the  case 
here  given  the  mirror  surface  of  the  Bil3  single  crystals  investigated  here.  Therefore,  this 
simplified  solution  of  absorption  coefficient  will  not  truly  represent  the  band  gap 
characteristics  of  Bil3  as  key  features  may  be  lost  in  the  (ahv)1/2  versus  hv  plot  when  a 
approaches  zero. 

Nevertheless,  it  is  clearly  observed  that  in  the  UV-Vis  transmission  spectra  the 
absorption  edge  of  Bil3  red-shifts  with  increasing  temperature  (Figure  26).  The 
temperature  coefficient  between  298  and  330  K  was  calculated  based  on  the  plot  in 
Figure  27a  and  the  value  is  about  -1 .1 8  x  1 0'3  eV/K.  This  temperature  effect  has  also 
been  reported  by  other  researchers,  as  indicated  in  Table  4-1.  This  red-shifting  of  the 
band  gap  with  temperature  has  been  generally  found  in  many  semiconductor 
materials.62'64  An  empirical  fit  has  been  experimentally  determined  to  describe  the  band 
gap  as  a  function  of  temperature:65 


r/T2 

E,(T)=  E,(0)~—  (13) 

where  Eg(0)  is  the  band  gap  at  absolute  temperature,  a  and  (3  are  the  fitting  parameters. 
Fitting  the  data  obtained  from  Figure  27a  gives  Eg(0)  =  2.06  eV,  a  =  1 .36x1 0‘3  eV/K  and 
P  =  52  K.  The  band  gap  at  absolute  0  K  is  then  determined  to  be  2.06  eV,  which  is 
about  0.2  eV  smaller  than  the  theoretical  value.  If  aT/(T+(3)  is  taken  as  temperature 
coefficient  of  the  band  gap,  this  coefficient  increases  with  increasing  temperature. 
Different  temperature  coefficients  with  close  values  have  been  reported  for  Bil3.  For 
example,  Komatsu  and  Kaifu  obtained  temperature  coefficient  of -1.8  x  10'4  eV/K 
between  6  and  77  K,  and  -3.4  x  10'4  eV/K  between  77  and  286  K.  Vashchenko  and 
Timofeev  reported  a  coefficient  of  -9.2  x  10'4  eV/K  between  85  and  293  K.  As  a 
comparison,  temperature  coefficients  of  Bil3  were  calculated  using  above  fitted 
parameters  and  listed  in  Table  7. 


Table  7.  Temperature  coefficient  of  Bil3  calculated  based  on  the  fitting  parameters. 


Temperature  (K)  6 

77 

85 

286  293 

Temperature 

-8.1  xIO4 

-8.4  x  104 

-1.2  x  KV5  -1.2  X1CT3 

Coefficient  -I.4xl04 

(eV/K) 

It  is  shown  that  temperature  coefficients  obtained  in  this  study  have  slightly  higher 
values,  yet  the  range  of  the  values  have  good  agreement  with  those  reported 
previously. 
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3.2.10  Effects  of  Impurities  on  Band  Gap 


Band  gap  narrowing  is  an  effect  commonly  observed  in  heavily  doped  semiconductor 
materials.66"68  The  physical  mechanisms  for  band  gap  narrowing  fall  into  two  categories: 
dopant  atoms  interacting  with  the  semiconductor  and  modifying  its  band  structure,  and 
strain  introduced  by  the  doping  process.69 The  band  gap  narrowing  effect  has  a  general 
form  of 


AEg  »  AN^  ,  (14) 

where  A  is  constant,  and  N  is  the  dopant  concentration.69  For  Si,  the  band  gap 
narrowing  is  about  0.02  eV  with  1017  cm'3  dopant  concentration.70  If  the  impurities  in  Bil3 
single  crystal  are  considered  as  an  extrinsic  dopant  and  the  concentration  is  high,  a 
similar  band  gap  narrowing  effect  is  expected. 

The  impurity  levels  of  the  Bil3  powder  along  with  the  crystals  grown  in  Pyrex  glass  and 
quartz  ampoules  are  reported  in  Table  8.  The  concentration  was  calculated  based  on 
the  mole  ratio  between  the  impurity  element  and  Bil3.  The  concentration  of  Pb  in  the 
powder  and  crystals  was  below  the  detection  limits  of  ICP-AES. 

Table  8.  Impurity  Concentrations  in  Bil3  Powder  and  Crystals. 


Aq 

Cu 

Fe 

Ni 

Pb 

Si 

Raw 

powder 

3.4±0.60 

0.8±0.02 

7.1±1.86 

0.1±0.01 

<1 

7.2±0.71 

Impurity  concentrations  of  crystal  grown  in  Pyrex  glass  amDoule 

Distance 
from  the 
tip  (mm) 

Ag 

CU 

Fe 

Ni 

Pb 

Si 

0 

2.2±0.57 

0.3±0.01 

5.5±0.05 

0.2±0.01 

<1 

4.8±0.38 

5 

4.6±0.55 

0.3±0.01 

6.1±0.08 

0.1±0.01 

<1 

6.5±0.46 

10 

2.4±0.10 

0.5±0.01 

6.0±0.24 

0.2±0.01 

<1 

5.9±0.78 

15 

2.9±0.13 

0  2±0.01 

6.6±0.12 

0.1±0.01 

<1 

6.3±0.53 

20 

1.3±0.03 

0.2±0.01 

4.3±0  18 

0.1±0.01 

<1 

6.5±0.38 

25 

2.6±0.06 

0.1  ±0.01 

5.3±0.08 

0.1±0.01 

<1 

4.6±0.34 

Impurity  concentrations  of  crystal  grown  in  guartz  ampoule. 

0 

2.4±0.1 1 

0.2±0.01 

5.6±0.21 

0.3±0.01 

<1 

4.2±0.15 

5 

2.0±0.17 

0.2±0.01 

4.2±0.1 1 

0.3±0.02 

<1 

3.5±0.23 

10 

3.1  ±0.08 

0.2±0.01 

4.9±0.15 

0.1±0.01 

<1 

4.7±0.36 

It  was  observed  that  the  impurity  concentrations  of  the  grown  crystal  maintained  the 
same  level  as  Bil3  powder.  The  crystals  grown  in  Pyrex  and  quartz  ampoules  had 
similar  impurity  concentrations,  which  suggested  that  crystal  growth  in  Pyrex  glass 
ampoule  was  sufficient  to  maintain  the  impurity  level  of  the  starting  material.  This  is  a 
great  processing  advantage  given  the  easier  ampoule  vacuum  sealing  of  Pyrex  when 
compared  to  quartz. 
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The  sum  of  the  impurity  concentrations  for  all  the  analytes  were  around  20  ppm  in  the 
crystal.  Given  that  Bil3  has  an  atomic  concentration  of  2.37  x  1022  cm'3,  20  ppm  impurity 
concentration  was  converted  into  an  atomic  concentration  of  4.74  x  1017  cm-3.  This  was 
considerably  high  doping  concentration  if  all  the  impurity  elements  were  “doped”  into  the 
crystal.  It  was  therefore  reasonable  to  expect  that  the  band  gap  of  Bil3  crystal  was 
narrowed  by  the  impurities. 

The  impurity  distribution  in  the  crystal  grown  in  Pyrex  glass  ampoule  was  uniform  along 
the  axis  of  the  crystal  ingot,  as  shown  in  Figure  28,  which  suggested  that  the  impurities 
was  hardly  segregated  during  MVB  crystal  growth. 


Figure  28.  Impurity  distribution  along  the  Bil3  crystal  grown  in  Pyrex  ampoule. 

The  band  gap  measured  at  room  temperature  showed  close  value  to  most  of  the 
literature  listed  in  Table  6  except  for  reference  39  and  71 .  Regarding  to  reference  71 , 
one  should  be  noted  that  the  values  were  actually  taken  from  the  peak  energy  of  the 
fundamental  absorption  edge,  which  should  be  higher  than  the  intrinsic  band  gap  as  the 
value  is  taken  at  the  maximum  absorption  coefficient  a  instead  of  at  a  approaching  to 
zero.  In  other  words,  the  values  reported  in  reference  71  are  related  to  band  gap  but 
with  higher  energy.  They  are  not  the  true  band  gap  values,  but  their  change  with 
temperature  represents  the  temperature  dependence  feature  of  band  gap. 

The  band  gap  in  reference  39  was  measured  by  ellipsometry.  A  comparison  of 
spectroscopic  ellipsometry  and  transmission  and  reflection  spectroscopy  measurements 
indicates  the  relative  strengths  and  weaknesses  of  the  two  techniques.  Spectroscopic 
ellipsometry  has  proven  to  be  useful  in  extracting  the  material  microstructure  in  the  form 
of  surface  features,  the  index  of  refraction  below  the  band  gap,  and  the  critical  point 
features  above  the  direct  gap.  This  technique,  however,  lacks  the  sensitivity  when  the 
absorption  coefficient  is  lower  than  a  ~  104  cm"1.  When  the  absorption  coefficient  is 
extracted  using  transmission  and  reflection  spectroscopy,  sensitivity  is  retained  to  a  ~ 
102  cm'1.  Transmission  and  reflection  measurements  generally  lack  sensitivity  to  the 
material  microstructure,  which  can  lead  to  erroneous  assignment  of  the  above-gap 
critical  point  transitions.  A  combination  of  these  two  techniques  allows  for  the  accurate 
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determination  of  the  material  microstructure  and  of  the  optical  properties  across  the 
given  spectral  range.  From  reference  72,  it  is  clearly  apparent  that  when  the  band  gap 
characteristics  are  above  a  ~  104  cm-1  either  technique  is  useful  for  extraction  of  the 
direct  or  indirect  gap,  but  when  the  characteristics  are  lower  than  a  ~  104  cm'1 
transmission  and  reflection  measurements  should  be  used.  Thus  the  different  band  gap 
values  obtained  by  UV-Vis  and  ellipsometry  can  be  attributed  to  the  different  range  of 
absorption  coefficient  used  for  the  data  fitting,  which  is  determined  by  the  nature  of  the 
techniques.  In  the  case  of  Bil3,  UV-Vis  is  more  accurate  to  determine  the  band  gap  as 
its  band  characteristics  are  lower  than  a  ~  104  cm'1. 

The  decrease  of  band  gap  value  due  to  increasing  temperature  and  high  impurity 
concentration  has  been  observed  in  Bil3.  However,  the  effects  are  rather  small,  such 
that  they  cannot  match  the  difference  between  experimental  and  theoretical  values. 
Theoretical  DFT  calculation  gave  an  indirect  band  gap  value  of  2.28  eV.72  The 
experimentally  determined  value  was  2.06  eV  at  zero  temperature  by  fitting  the  data 
obtained  from  298  to  330  K.  The  band  gap  narrowing  caused  by  impurities  has  been 
theoretically  studied  using  DFT  calculation,  which  gave  a  0.3  eV  decrease  of  band  gap 
for  Cu  impurity  with  concentration  of  ~1020  cm'3.72  The  theoretical  study  also  showed 
that  not  every  impurity  element  will  cause  the  band  gap  narrowing  of  Bil3.  For  instance, 
Si  was  observed  to  have  no  effect  on  the  band  structure  modification.  The  band  gap 
narrowing  in  Bil3  should  be  far  less  than  0.3  eV,  since  the  impurity  concentration  in  the 
crystal  (~1017  cm'3)  is  much  smaller  than  the  value  used  in  DFT  calculation. 
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4. 


Detector  Fabrication  and  Testing 


4.1  Detector  Fabriaction 

For  detector  fabrication,  a  metal  electrode  has  to  be  deposited  on  the  surface  of  Bib. 

The  proper  choice  of  the  electrode  materials  is  extremely  important  for  the  success  of 
the  detector.  Au  and  Pd  have  both  been  used  as  electrode  materials  for  Bil3  detectors.1' 

3  However,  the  properties  of  the  electrical  contacts  have  never  been  reported.  It  is 
important  to  understand  how  the  metal  contact  interacts  with  the  substrate.  The 
reactivity  of  the  metal-semiconductor  interface  should  be  used  as  a  guide  for  the 
selection  of  contact  material.  X-ray  photoelectron  spectroscopy  (XPS)  has  been  widely 
used  to  explore  the  metal-semiconductor  interfaces  with  great  success  at  explaining 
and  predicting  the  reactivity  of  the  interfaces.73'74  In  this  study,  different  electrode 
materials  of  Au,  Pd,  and  Pt  have  been  sputter  coated  on  the  surfaces  of  Bil3  single 
crystals.  The  metal-Bil3  interface  has  been  studied  by  XPS  and  AFM.  The  resistivity  of 
Bil3  single  crystals  with  different  contact  materials  have  been  measured  and  compared. 

4.1.1  XPS 

4.1. 1.1  Principles  of  XPS 

XPS  is  an  electron  spectroscopic  method  that  uses  monoenergetic  soft  X-rays  to  eject 
electrons  from  the  inner  shell  orbitals  of  a  sample  and  analyzes  the  energy  of  the 
detected  electrons.75  In  this  study,  the  X-ray  source  is  a  monoenergetic  Al  Kc  (1486.7 
eV)  X-ray.  The  kinetic  energy  of  an  ejected  electron  is  given  by 

Ek=hv-Eb-<ps  (15) 

where  hv  is  the  energy  of  the  photon,  Eb  is  the  electron  binding  energy  which  is  the 
energy  difference  between  the  Fermi  level  and  the  energy  of  the  atomic  orbital  from 
which  the  electron  originates,  and  <ps  is  the  work  function  of  the  spectrometer.  Since  an 
atom  has  multiple  orbitals,  the  corresponding  photoelectron  spectrum  has  a  variety  of 
peaks  that  are  named  according  to  the  orbital  (I  =  0,  1 , 2,  3... denoted  as  s,  p,  d,  f...)  and 
spin  (s  =  1/2)  quantum  numbers.  The  total  momentum  of  the  photoelectrons  (J  =  I  +  s)  is 
included  in  the  labeling  of  a  measured  photoelectron  peak  (e.g.,  Bi  4f7/2  where  I  +  s  =  3 
+  1/2  =  7/2).  The  electron  binding  energy  of  an  element  varies  with  different  compounds 
due  to  the  change  of  electronic  environment.  Therefore  XPS  can  be  used  to  identify  the 
chemical  state  of  the  materials  being  analyzed.  Since  the  energy  of  the  X-ray  source  is 
low,  the  penetration  depth  of  the  photons  is  on  the  order  of  micrometers.  Only  those 
electrons  that  originate  within  tens  of  angstroms  below  the  solid  surface  can  escape  the 
surface  without  appreciable  energy  loss,  which  makes  XPS  majorly  used  as  a  surface 
analytical  technique. 

4.1. 1.2  Sample  Preparation  for  XPS 
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The  (001)  surface  of  a  Bil3  single  crystal  was  cleaved  in  air  by  peeling  off  the  top  layer 
of  the  crystal  using  Scotch  tape.  The  cleaved  surface  was  sputtered  with  Au,  Pd,  and  Pt 
metals  with  thickness  of  around  1  nm  using  Kurt  J.  Lesker  CMS-18  multi  target  sputter 
system.  The  sputtering  was  conducted  in  a  vacuum  chamber  with  a  pressure  of  10'8  tori 
at  a  deposition  speed  of  around  0.5  nm/s  at  room  temperature.  A  freshly  cleaved 
surface  of  a  Bil3  single  crystal  was  used  as  a  reference.  X-ray  photoelectron  spectra 
were  collected  by  an  Omicron  XPS  system  with  an  Al  Ka  monochromatic  X-ray  (1486.7 
eV).  The  spectra  were  recorded  at  a  sampling  area  of  1 .3  mm2  and  a  take-off  angle  of 
55°  with  pass  energy  of  22  eV  and  step  size  0.05  eV.  The  chamber  pressure  was  kept 
below  5  x  10'10  mbar  during  measurement.  The  instrument  resolution  was  <  0.1  eV.  The 
peak  energies  were  calibrated  using  a  Ag  3d5/2  peak  at  368  eV.  The  studied  elements 
included  bismuth,  iodine,  gold,  palladium,  platinum,  carbon,  and  oxygen. 

4.1.2Atomic  Force  Microscopy  (AFM) 

To  study  the  surface  morphology  of  the  metal  contacts,  the  sample  surfaces  were 
examined  by  AFM  with  an  Asylum  MFP-3DTM  Research  AFM.  The  AFM  tip  from  Veeco 
has  a  triangular  cantilever  and  a  nominal  tip  radius  of  20  nm.  The  normal  force  constant 
of  the  tip  cantilever  is  0.58  N/m.  The  topography  images  were  acquired  in  contact  mode 
at  a  fixed  applied  load  of  several  nN. 

4.1.3  l-V  Characterization 

The  resistivity  of  Bil3  single  crystals  with  different  electrodes  was  studied  by  l-V 
characterization  using  an  Agilent  4155B  semiconductor  parameter  analyzer.  Samples 
were  prepared  by  sputtering  Au,  Pd,  and  Pt  on  the  freshly  cleaved  surface  resulting  in 
an  electrode  thickness  of  ~200  nm.  The  sputtering  was  conducted  at  the  same 
conditions  as  described  for  the  XPS  samples.  Pt  wires  were  attached  to  the  electrodes 
using  silver  paste  (SPI®  Supplies)  for  electrical  connection. 

4.1.4  Results  and  Discussion 
4.1. 4.1  The  XPS  Spectra 

The  XPS  core-level  spectra  of  all  the  electrode  elements  detected  are  shown  in  Figure 
29.  The  binding  energies  of  each  element  in  different  compounds  are  given  in  Table  9. 
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Figure  29.  XPS  spectra  of  Bi  4f  and  I  3d  for  the  freshly  cleaved  Bil3  (001)  surface  and  surface  with  Au,  Pd, 

and  Pt  electrodes. 

Table  9.  Observed  binding  energies  of  the  several  elements  observed  within  the  XPS  spectra. 


Peak 

Bib  sample 

Binding  energy  (eV) 

Au  samole  Pd  samtfe 

Pt  samDle 

Bi  4f7/2 

159.2 

159.1 

158.6  157.1 

150.6  157.3 

I3CJ5/2 

619.4 

619.1 

619.3 

619.0 

C  is 

284.7 

284.6 

284.0 

284.3 

O  is 

532.2  529.9 

533.0 

532.1 

Au  4f7/2 

04.0 

Pd  3d5/2 

335.0  336.5 

Pt  4f7/2 

71.1 

In  the  following  discussion,  the  sample  without  deposited  metal  is  called  Bil3  sample, 
and  the  samples  with  deposited  metallic  contacts  are  named  Au,  Pd,  and  Pt  samples 
respectively.  A  quick  survey  of  the  XPS  spectra  reveals  that  the  Au  sample  maintains 
the  chemical  environment  of  Bil3,  while  Pd  and  Pt  show  interactions  with  the  substrate 
The  detailed  discussion  for  each  sample  is  given  in  the  following  sections. 
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4.1.4.2  Au  Sample 


The  XPS  spectra  of  the  Au  sample  show  peak  positions  of  Bi  4f  and  I  3d  very  close  to 
those  of  the  Bil3  sample;  and  the  Au  4f  spectrum  exhibits  the  same  peak  positions  with 
that  of  a  pure  gold.  The  results  suggest  that  Bi,  I,  and  Au  maintain  their  respective 
chemical  state  and  there  is  no  chemical  interaction  between  Au  and  Bil3.  The  C  Is  and 
O  Is  spectra  are  shown  in  Figure  30. 


Figure  30.  Deconvolution  of  XPS  spectra  for  the  AU  sample  peak  (a)  C  1  s  and  (b)  O  1  s. 

The  C  Is  spectrum  of  the  Au  sample  can  be  deconvolved  into  two  peaks  with  the  peak 
at  284.5  eV  assigned  to  amorphous  carbon76  and  the  peak  at  285.7  eV  assigned  to 
carbon  species  bound  to  oxygen77  79(Figure  30a).  Two  distinct  peaks  are  observed  in  O 
Is  spectrum  (Figure  30b).  The  peak  at  532.2  eV  is  attributed  to  oxygen  atoms  bound  to 
carbon.77  The  peak  at  529.9  eV  could  be  assigned  to  oxygen  atoms  in  Au203  (530.2 
eV)58  or  atomic  O  adsorbed  at  Au  surface  (530.2  eV).77  Since  there  is  no  indication  of 
any  Au-0  bond  in  the  Au  4f  spectrum,  this  peak  is  assigned  to  atomic  O  adsorbed  at  the 
Au  surface. 

4.1. 4.3  Pd  Sample 

Peak  shifting  of  Bi  4f  and  Pd  3d  is  observed,  which  indicates  the  existence  of  an 
interfacial  reactivity.  The  iodine  peaks  remain  the  same  as  the  Bil3  sample.  The 
deconvolution  of  Bi  4f,  Pd  3d,  C  Is,  and  Ols  spectra  are  shown  in  Figure  31.  The 
properties  of  the  fitted  peaks  and  their  assignments  are  summarized  in  Table  10. 
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Figure  30.  Deconvolution  of  XPS  spectra  of  the  Pd  sample  for  (a)  Bi  4f,  (b)  Pd  3d,  (c)  C  Is  and  (d)  O  Is. 


Table  10.  Peak  Properties  of  Bi  4f,  Pd  3d,  Pd  3p,  C  Is,  and  O  Is  Observed  in  the  Pd  Sample. 


Spectrum 

Position 

(eV) 

Element 

FWHM 

(eV) 

Area 

ASF 

Bi  4f 

157.1 

Bi  4f7/2 

0.87 

716 

7.4 

Bi4f 

158.5 

Bi  4f7/2 

0.87 

460 

7.4 

Bi  4f 

159.2 

Bi  4f7/2 

0.87 

152 

7.4 

Bi  4f 

162.4 

Bi  4f5/2 

0.87 

537 

7.4 

Bi  4f 

163.8 

Bi  4f5/2 

0.87 

345 

74 

Bi  4f 

164.4 

Bi  4f5/2 

0.87 

114 

74 

Pd  3d 

335.0 

P  d  3d  5/2 

1.03 

4512 

46 

Pd  3d 

336.5 

Pd  3d  5/2 

1.03 

6796 

4.6 

Pd  3d 

340.3 

Pd  3d  3/2 

1.03 

2978 

4.6 

Pd  3d 

341.7 

P  d  3d  3/2 

1.03 

4486 

4.6 

C  IS 

283.9 

C  IS 

1.38 

871 

0.25 

C  is 

285.6 

CIS 

1.38 

143 

0.25 

O  IS 

533.1 

Pd  3p3/2 

3.G9 

5862 

143 

OIS 

5321 

OIS 

1.80 

378 

0.66 

Corrected 

Area 

97 

62 

21 

73 

47 

15 

961 

1477 

647 

975 

3464 

572 
4099 

573 


Assignment 

Bi 

Bil3 

Bi 

BiU-x 

Bib 

Pd 

Pdb 

Pd 

Pdl2 

amorphous 
carbon 
C  bound  to 
O 
Pd 

O  bound  to 
C 


In  the  Bi  4f7/2  region,  three  peaks  are  identified  at  157.1,  158.5,  and  159.2  eV,  which 
are  consistent  with  the  literature  values  assigned  to  bismuth  atoms  in  metallic  bismuth,80 
Bi203,81  and  Bil375  respectively.  However,  the  assignment  of  Bi203  is  not  plausible  as 
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the  corresponding  amount  of  oxygen  is  not  observed  in  the  O  Is  spectrum.  The  peak  at 
158.5  eV  is  possibly  due  to  bismuth  atoms  in  a  reduced  iodine  environment  (Bil3.x  with  x 
<  3).  The  Pd  3d5/2  peak  deconvolutes  into  two  peaks  at  335.0  and  336.5  eV,  which  are 
assigned  to  palladium  atoms  in  metallic  palladium75  and  Pdl22  respectively. 

The  C  spectrum  of  the  Pd  sample  has  two  peaks  with  the  one  at  283.9  eV  assigned  to 
amorphous  carbon  and  the  one  at  285.6  eV  assigned  to  carbon  species  bound  to 
oxygen.75  The  O  Is  spectrum  is  dominated  by  Pd  3p3/2peak  at  533.1  eV.75  The  intensity 
of  this  peak  (area  under  the  peak)  is  consistent  with  the  peak  intensity  in  Pd  3d 
spectrum,  which  validates  this  peak  assignment.  The  identification  of  oxygen  peak  in  O 
Is  spectrum  is  difficult  as  the  majority  peak  intensity  is  contributed  by  Pd  atoms.  The 
peak  is  fitted  at  532.1  eV  with  a  fixed  1 :1  intensity  ratio  to  the  corresponding  C  Is  peak 
at  285.6  eV  in  C  Is  spectrum. 

XPS  analysis  indicates  the  existence  of  metallic  Bi  and  Pdl2  at  the  Pd-Bil3  interface, 
which  suggests  that  Pd  reacts  with  Bil3.  The  majority  of  the  bismuth  elements  on  the 
surface  exist  in  the  form  of  metallic  Bi  (~54%)  and  Bii^x  (~35%),  with  a  small  portion 
(~11%)  that  remains  as  Bil3.  The  reactivity  between  Pd  and  iodine-rich  crystals  has 
been  previously  reported.82  However,  Pd  is  still  widely  used  as  the  electrode  contacts 
for  iodide  compound  materials  such  as  Hg!2  and  Pbl2.22,83-84  There  is  no  discussion  on 
how  this  Palladium-iodine  reactivity  will  affect  the  detector  performance.  This  suggests 
that  a  careful  re-evaluation  of  the  electrode  contacts  for  these  materials  is  necessary. 

4.1. 4.4  Pt  Sample 

Similar  to  the  Pd  sample,  peak  shifting  is  also  observed  in  the  Bi  4f  spectrum  in  the  Pt 
sample.  The  Pt  4f  peaks  exhibit  an  asymmetric  feature,  indicative  multiple  chemical 
states  for  the  Pt  atoms  found  in  the  interfacial  region.  However,  this  asymmetry  is  not 
due  to  reactivity  therefore  no  attempt  is  made  to  deconvolve  the  peaks.  The  similar 
peak  feature  has  been  previously  observed  in  Pt  nanoparticles  on  pyrolytic  graphite 
surfaces.85  Since  Pt  was  reported  to  have  different  electronic  configurations  for  surface 
and  bulk  atoms,86  the  asymmetric  peak  was  explained  due  to  the  surface  Pt  atoms  in 
the  Sd^s1  configuration,  bulk  Pt  atoms  in  the  5d86s2  configuration,  and  surface  Pt 
absorbed  oxygen.85  The  peak  properties  and  their  assignments  are  summarized  in 
Table  11. 
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Table  1 1 .  Peak  properties  of  Bi  4f,  Pd  3d,  Pd  3p,  C  Is,  and  O  Is  observed  in  the  Pt  sample. 


Spectrum 

Position 

(eV) 

Element 

FWHM 

<eV) 

Area 

ASF 

Corrected 

Area 

Assignment 

Bi  4f 

157.1 

Bi  4f7/2 

0.86 

2303 

74 

311 

Bi 

Bi  4f 

158.3 

Bi  4f7/2 

0.86 

1501 

7.4 

203 

Bib-* 

Bi  4f 

158.9 

Bi  4f7/2 

0.86 

891 

7.4 

120 

Bil3 

Bi  4f 

162.4 

Bi  4f5/2 

0.86 

1727 

7.4 

233 

Bi 

Bi  4f 

163.6 

Bi  4f5/2 

0.86 

1125 

7.4 

152 

Bilj-x 

Bi  4f 

164.2 

Bi  4f5/2 

0.86 

669 

7.4 

90 

Bib 

13d 

618.6 

1  3d5/2 

1.12 

10385 

6 

1731 

Bib 

13d 

619.4 

1  3d5/2 

1.12 

4365 

6 

728 

b 

13d 

630.1 

1  3d3/2 

1.12 

7555 

N/A 

Bib 

13d 

630.9 

1  3d3'2 

1.12 

2550 

N/A 

b 

C  is 

284.0 

C  IS 

1.46 

884 

0.25 

0.25 

3536 

1036 

amorphous 

carbon 

C  bound  to 

C  Is 

285.9 

C  IS 

1.46 

259 

O 

0.66 

280 

O  bound  to 

Ols 

529.5 

O  IS 

1.59 

185 

Pt 

0.66 

945 

O  bound  to 

CMS 

532.0 

O  IS 

1.59 

624 

C 

The  assignment  for  Bi  4f  spectrum  is  similar  to  that  in  Pd  sample.  The  Bi  4f7/2  peaks  at 
157.1,  158.3,  and  158.9  eV  are  assigned  to  bismuth  atoms  in  metallic  bismuth,  Bil3.x, 
and  Bil3  respectively.  The  existence  of  l2  is  evidenced  by  the  peak  at  619.4  eV  in  I  3d 
spectrum.75  The  two  peaks  in  C  Is  spectrum  are  assigned  to  amorphous  carbon  at 
284.0  eV  and  carbon  species  bound  to  oxygen  at  285.9  eV.  The  0  1s  spectrum 
deconvolves  into  two  peaks  with  the  one  at  529.5  eV  assigned  to  oxygen  bound  to  Pt 
and  the  one  at  532.0  eV  assigned  to  oxygen  bound  to  carbon.  The  uphill  tail  in  the  low 
energy  side  of  O  Is  spectrum  is  the  trail  of  Pt  4p3/2  peak  at  520  eV.75 

The  XPS  results  indicate  that  Pt  is  not  good  as  the  contact  material  for  Bil3  either.  The 
effect  of  Pt  on  the  surface  of  Bil3  is  different  from  that  of  Pd.  Here  the  presence  of  Pt 
causes  the  reduction  of  Bil3.  But  unlike  Pd,  Pt  is  not  involved  in  the  reaction.  The 
function  of  Pt  in  this  case  is  similar  to  that  of  a  catalyst. 

4.1.5  AFM  Analysis 

AFM  characterization  has  been  used  to  check  the  distribution  of  the  sputtered  metals  on 
the  surface  of  Bil3  (Figure  31).  It  is  observed  that  Au  fully  covers  the  Bil3  surface  while 
Pd  and  Pt  form  islands,  which  indicates  a  completely  wetting  of  Au  on  the  Bil3  surface 
while  non-wetting  characteristic  of  Pd  and  Pt.  The  islands  on  the  Pt  sample  surface  are 
more  sparsely  distributed  with  higher  height  compared  to  the  Pd  sample.  To  verify  if  the 
islands  are  formed  by  the  deposited  metals  or  by  other  species  produced  in  the  metal- 
Bil3  interaction,  the  volume  of  the  islands  is  calculated  using  an  image  processing 
software  (Scanning  Probe  Image  Processing,  SPIPTM)  by  analyzing  the  scanned 
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histogram  of  the  surface.  The  calculated  island  volume  for  Pd  sample  is  4.6  x  106  nm3 
on  a  4  pm2  area,  which  is  corresponding  to  a  thickness  of  1.15  nm  if  the  Pd  island 
volumes  were  evenly  spread  on  the  surface.  This  matches  well  with  the  experimentally 
deposited  thickness  of  1  nm  and  suggests  that  the  islands  are  mainly  formed  by  Pd. 
The  island  volume  on  Pt  sample  is  calculated  to  be  7.7  x  108  nm3  on  a  100  pm2  area, 
which  is  corresponding  to  a  deposited  thickness  of  7.7  nm.  This  result  has  a  large 
deviation  from  the  experiment  (1  nm),  which  suggests  that  species  generated  by  the 
reaction,  i.e.  I2,  may  deposit  on  the  surface. 


(c) 


Figure  31.  AFM  image  and  3D  view  of  topography  for  the  (a)  Au  sample,  (b)  Pd  sample,  and  (c)  Pt 

sample. 

4.1.6  l-V  Characterization 

The  resistivity  of  Bil3  single  crystals  has  been  measured  with  Au,  Pd,  and  Pt  electrodes. 
The  crystal  samples  were  cut  from  adjacent  positions  of  a  large  single  crystal  to 
minimize  the  sample  variation.  The  l-V  curves  for  crystals  with  different  electrodes  are 
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shown  in  Figure  32.  It  is  observed  that  the  crystal  with  Au  electrodes  has  the  highest 
resistivity,  and  the  one  with  Pd  electrodes  shows  the  lowest,  as  summarized  in  Table 
12. 


Figure  32.  Typical  l-V  characteristic  curves  for  Bil3. 

Table  12.  Sample  Information  and  Resistivity  of  Bil3  Single  Crystals  for  l-V  Characterization 


Sample  # 

Electrode 

Surface^  Area 

(cm2) 

Thickness  (cm) 

Resistivity 

10s  (Q-cm) 

1 

Au 

0.260 

0.075 

1.85 

2 

Pd 

0.167 

0.075 

0.28 

3 

Pt 

0.251 

0.076 

0.61 

The  similar  results  were  also  reported  by  Fornaro  et  at.  on  Bil3  crystals  grown  by  PVT 
method.3  Given  that  Bi3+  ions  are  reduced  to  metallic  Bi  at  the  Pd-Bil3  interface,  the  Bi 
atoms  could  drift  across  the  crystal  body  under  an  external  electric  field.  This  will 
decrease  the  effective  thickness  between  the  electrodes.  Consequently,  the  measured 
resistivity  is  smaller  than  the  actual  value. 

Based  on  the  XPS  results,  Au  is  expected  to  be  physically  contacted  with  Bil3,  the 
measured  resistivity  thus  represents  a  more  accurate  value  of  this  material.  The  sample 
with  Pt  electrodes  shows  resistivity  in  between,  which  is  possibly  due  to  the  drift  of  Bi  + 
ions  resulting  from  less  binding  of  iodine. 

4.2  Radiation  Sensing 

We  cut  Bil3  crystals  using  a  diamond  wire  saw  purchased  for  this  project  from  MTI 
Corporation.  This  wire  saw  has  variable  pressure  and  saw  rotation  speed  which  allows 
us  to  investigate  for  cutting  and  surface  damage  due  to  the  detector  fabrication  process 
An  example  of  one  of  our  first  cut  detector  crystals  is  shown  in  Figure  33.  Detectors  of 
various  with  Pd  electrodes  were  fabricated  for  testing. 
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Figure  33.  Left:  Cut  Bil3  crystal  prepared  for  electrode  (Pd)  deposition.  Right:  Fabricated  detector  after 

electrode  deposition  and  wire  bonding. 

To  prepare  for  testing,  we  built  two  test  boxes  for  the  detectors.  We  measured  signals 
from  the  anode  and  cathode  electrodes.  Amptek  A250  charge  sensitive  preamplifiers 
were  installed  into  a  light-tight  aluminum  box.  We  also  developed  a  data  acquisition 
system  to  be  used  in  detector  testing.  This  program,  built  using  LabView,  has  the 
capability  to  digitize  the  preamplifier  signals  and  perform  acquire  both  spectral  and 
timing  information  from  the  detectors.  A  screen  shot  of  the  data  acquisition  interface  is 
shown  in  Figure  34.  The  system  was  tested  with  LaBr3  crystals  in  our  laboratory. 
Unfortunately,  the  sampling  speed  of  the  data  acquisition  card  is  1  MHz  per  channel  for 
4  channels.  This  sampling  speed  may  not  be  able  to  accurately  measure  from  the 
preamplifier.  Data  from  shaping  amplifiers  showed  much  improved  energy  accuracy  and 
improved  energy  resolution.  Thus,  we  modified  the  program  to  accept  both  signals  from 
the  preamplifier  and  shaping  amplifier  to  ensure  that  we  can  collect  all  information  from 
the  detector. 


Figure  34.  Screenshot  showing  front  display  for  LabView  data  acquisition  program  to  collect  spectral  and 
timing  information  from  detectors  (Left).  Test  box  with  preamplifiers  and  detector  (Right). 
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4.2.1  Measurement  of  Electron  Mobility 


The  mobility  of  electrons  and  holes  is  one  property  that  helps  describe  the  charge 
transport  characteristics  in  a  semiconductor  detector.  A  detector  with  high  and  uniform 
charge  carrier  properties  throughout  the  material  will  generally  perform  better  as  a 
radiation  detector  and  spectrometer.  In  planar  detectors  with  ohmic  contacts,  the 
electric  field  should  be  relatively  uniform.  The  mobility  (p)  of  electrons  and  holes  can  be 
estimated  by  measuring  the  drift  time  (W)  from  an  interaction  in  a  detector.  To  estimate 
the  mobility  of  electrons,  one  can  use  alpha  particles  that  will  be  stopped  in  a  few  pm, 
irradiated  from  the  cathode  side  of  the  detector,  forcing  electrons  to  drift  the  entire 
thickness  (D)  of  the  detector.  In  this  case,  the  mobility  of  electrons  can  be  calculated  by 


Et drift, e 


where  E  is  the  strength  of  the  electric  field  (V/cm),  and  the  subscript  e  represents 
electrons.  We  collected  a  pulse  waveform  from  a  Bib  detector  (Figure  35).  By 
measuring  the  rise  time  of  the  signal,  we  get  the  drift  time  of  the  electrons. 
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Figure  35.  Pulse  waveform  from  an  Am-241  alpha-particle  irradiated  from  the  cathode.  By  measuring  the 
drift  time  of  the  electrons  creating  the  signal,  a  mobility  for  electrons  can  be  estimated. 

The  drift  time  for  the  electrons  in  the  above  waveform  was  estimated  to  be 
approximately  250  ns  (0.25  ps).  Thus,  using  the  detector  thickness  of  0.071  cm  and  the 
applied  voltage  of  70  V,  we  estimate  the  mobility  of  electrons  in  this  detector  to  be 
approximately  250±50  cm2A/-s.  The  uncertainty  in  this  value  is  based  on  several 
factors,  including  the  range  of  alpha  particles  in  Bib,  contributions  from  electronic  noise, 
and  uncertainty  in  drift  time.  It  should  also  be  noted  that  this  value  is  approximately  a 
factor  of  2  less  than  previously  reported  values  for  Bib  electron  mobility.  We  are 
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currently  investigating  reasons  for  the  discrepancy.  It  should  also  be  noted  that  the 
signal  strength  of  only  ~1  mV  was  much  lower  than  expected.  We  estimated  that  the 
signal  amplitude  should  have  been  about  40  mV.  Possible  reasons  for  the  poor  signal 
(and  additional  uncertainty  in  mobility)  include  dead  layer  formation  near  the  electrodes 
where  alpha  particle  would  be  stopped,  and  poor  electrode  contacts  where  the  signal 
size  would  be  reduced  due  to  weighting  potential  effects. 

4.2.2  X-Ray/Gamma-Ray  Detection 

The  detectors  were  tested  with  a  low-current  (0.1  mA)  60  kVp  x-ray  source.  The 
detectors  were  connected  to  an  Amptek  A250  preamplifier  and  were  connected  to 
standard  Ortec  NIM  electronics.  An  example  spectra  is  shown  on  the  right  in  Figure  36. 
Signals  from  the  preamp  (analyzed  on  an  oscilloscope)  were  almost  an  order  of 
magnitude  lower  than  expected  (<  0.5  mV  instead  of  8-10  mV).  This  small  signal 
suggests  that  improvements  in  surface  preparation  may  be  necessary  in  order  to 
increase  the  signal  from  the  detector. 


Figure  36.  60  kVp  X-ray  spectra.  A  1  cm2  collimated  beam  was  irradiated  on  the  cathode  side  of  the 

detector.  (-50  V  bias,  1  ps  shaping  time). 
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5. 


Conclusions 


The  increasing  demand  for  room  temperature  gamma-ray  detectors  encourages  the 
discovery  and  development  of  high  quality  compound  semiconductor  materials.  The 
required  material  properties  include  high  density  and  atomic  numbers  for  high  gamma- 
ray  detection  efficiency;  moderate  bandgap  between  1 .5  eV  and  3  eV  for  room 
temperature  operation  and  generation  of  large  number  of  electron-hole  pairs,  and  high 
resistivity  for  improved  signal  to  noise  ratio.  Bib  meets  all  these  criteria  thus  has  the 
potential  to  be  a  successful  material  for  room  temperature  gamma-ray  detection 
applications. 

Crystal  growth  is  the  most  challenging  and  time-consuming  stage  of  developing  a  new 
detector  material.  This  work  presented  the  efforts  of  growing  large  volume  Bil3  single 
crystals  using  MVB  method.  The  whole  growth  process  is  composed  of  different  stages 
including  heating  up  from  room  temperature  to  the  maximum  process  temperature, 
soaking  for  a  period  of  time,  cooling  down  to  growth  temperature,  crystal  growth,  and 
post-growth  cooling  down  to  room  temperature.  Each  stage  has  been  carefully  designed 
and  tested  to  ensure  the  smooth  growth  run.  The  crystal  growth  was  conducted  at  a 
total  of  nine  growth  conditions  (variation  of  growth  rate  and  temperature  gradient  at 
solid-liquid  interface),  among  which  only  three  conditions  produced  single  crystals.  The 
phase  purity  of  the  grown  single  crystals  was  verified  by  XRD.  The  quality  of  the  single 
crystals  was  characterized  by  EPD  measurement  and  X-ray  rocking  curve  w-scans.  For 
EPD  measurement,  the  etching  recipe  was  determined  to  be  10  sec  immersion  in  15% 
Kl  solution.  The  observed  EPD  values  for  all  the  crystals  were  on  the  order  of 
magnitude  of  ~104  cm'2,  which  was  comparable  to  other  compound  semiconductor 
crystals  such  as  Hgl2  and  lab  produced  CZT,  but  one  order  of  magnitude  higher  than 
the  commercial  CZT  crystals.  The  high  EPD  values  indicated  the  high  dislocation 
densities  in  the  crystals.  Among  the  crystals  grown  by  MVB  method,  the  one  grown  at 
10°C/cm  and  0.5  mm/h  showed  the  smallest  average  EPD  values.  The  rocking  curve  w- 
scan  peaks  for  all  the  crystals  showed  similar  FWHM  and  peak  characteristics.  Peak 
broadening  with  multiple  peaks  and  shoulders  were  observed.  The  poor  crystalline 
perfection  was  commonly  observed  among  the  iodide  compound  crystals  such  as  Hgl2 
and  Pbl2,  which  implied  that  the  imperfections  were  possibly  caused  by  the  residual 
strains  induced  during  the  cooling  of  the  crystals  due  to  the  softness  of  these  materials. 
In  addition,  the  deformation  during  handling  and  cutting  could  also  contribute  to  the 
broadening  of  the  rocking  curves.  Although  the  Bil3  single  crystals  grown  at  three 
different  conditions  did  not  show  distinct  differences  of  EPD  values  and  rocking  curve 
characteristics,  the  crystals  grown  at  10°C/cm  and  0.5  mm/h  were  selected  for  further 
studies  due  to  their  slightly  better  qualities.  A  single  crystal  with  a  dimension  of  18  x  13 
x  5  mm3,  the  largest  crystal  size  ever  reported  for  Bil3,  was  obtained  at  this  growth 
condition. 

The  cleaved  surface  of  Bil3  was  characterized  by  SEM  and  AFM  respectively.  The  SEM 
characterization  indicated  a  crack-free  and  macroscopically  smooth  surface.  The  AFM 
topographic  image  revealed  an  atomically  rough  (001)  surface  with  average  roughness 
of  0.486  nm.  The  irregular  edges  of  the  atomic  layers  were  attributed  to  the  weak  van 
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der  waal’s  forces  between  the  layers.  The  chemical  composition  at  the  surface  of  Bil3 
was  characterized  to  be  of  BH2.78  by  XPS.  The  iodine  deficiency  was  attributed  to  the 
dissociation  of  Bil3  into  Bi  and  I2  during  cooling  of  the  crystal,  and  the  evaporation  of 
iodine  when  exposed  in  air.  The  surface  oxidation  was  also  observed  for  Bil3  crystals 
exposed  in  air  for  over  two  weeks,  which  suggested  that  a  proper  sealing  or  inert  gas 
environment  is  necessary  for  the  crystal  storage. 

Through  UV-Vis  transmission  and  reflection  spectra,  the  band  gap  of  Bil3  was 
characterized  to  be  of  1 ,68  ±  0.09  eV  with  indirect  transition  characteristic.  The  band 
gap  value  of  Bil3  was  observed  to  decrease  with  increasing  temperature,  which  is 
commonly  observed  in  semiconductor  materials.  The  impurity  concentrations  in 
commercial  powder  and  Bil3  crystals  were  observed  to  be  as  high  as  1017  cm'3,  at  which 
the  band  gap  narrowing  effect  was  likely  to  happen.  The  band  gap  value  determined  in 
this  work  was  consistent  with  previous  studies  using  UV-Vis  spectra.  The  indirect  band 
transition  was  confirmed  by  the  most  recent  theoretical  calculation  on  the  band  structure 
of  Bil3  using  DFT.  The  large  difference  between  the  band  gaps  determined  by  UV-Vis 
spectra  and  ellipsometry  was  attributed  to  the  different  absorption  coefficient  a  used  to 
extract  the  band  gap.  For  Bil3,  band  gap  determined  by  UV-Vis  spectra  was  more 
accurate. 

Through  the  characterizations  on  Bil3  single  crystals  grown  in  this  work,  it  was  found 
that  the  crystals  have  high  impurity  concentrations  and  dislocation  densities.  The  crystal 
quality  is  not  comparable  to  the  commercial  CdTe  and  CZT  crystals.  While  the  material 
was  observed  to  be  sensitive  to  radiation  (alphas,  gamma-rays,  and  x-rays),  the  current 
material  was  not  stable  and  would  severely  polarize,  which  may  be  attributable  to  the 
high  impurity  concentrations  and  motion  of  impurities  at  high  voltage.  The  future  work 
should  focus  on  resolving  these  problems,  improving  the  purity  of  the  material,  and 
developing  Bil3  towards  a  room  temperature  gamma-ray  spectrometer. 
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6.  Additional  Information 

6.1  Papers  and  Presentations  Created  in  Whole  or  in  Part  by  This  Funding 

The  following  papers  and  presentations  were  made  possible  by  the  funding  received  by 

this  project  (in  whole  or  in  part).  In  addition  to  the  listed  publications,  at  least  one 

additional  paper  and  Ph.D.  thesis  associated  with  this  project  should  be  completed  in 

the  future. 

6.1.1  Publications 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Characterization  of  Bismuth 
Tri-Iodide  Single  Crystals  for  Wide  Band-Gap  Semiconductor  Radiation  Detectors,” 
accepted  in  Nuclear  Instruments  and  Methods  in  Physics  Research  Part  A,  201 1 . 

•  Qiu,  Wei,  “Growth  and  Characterization  of  Bismuth  Tri-Iodide  Single  Crystals  by 
Modified  Vertical  Bridgman  Method,”  Ph.D.  Thesis,  University  of  Florida,  2010. 

•  W.  Qiu,  A.  T.  Lintereur,  J.  C.  Nino,  and  J.  E.  Baciak,  “Growth  and  Characterization  of 
Bib  Single  Crystals  for  Room  Temperature  Gamma-Ray  Detection  Grown  by 
Modified  Bridgman  Method,”  submitted  to  Journal  of  Crystal  Growth,  2009. 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Growth  and  Characterization  of 
Bil3  Single  Crystals  for  Room  Temperature  Gamma-Ray  Detection  Grown  by 
Modified  Bridgman  Method,”  Hard  X-Ray,  Gamma-Ray  and  Neutron  Detectors 
Physics,  Proceedings  of  SPIE2,  Vol.  7449,  pp.  744958  1-8,  2009. 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Iodine  Based  Compound 
Semiconductors  for  Room  Temperature  Gamma-Ray  Spectrscopy,”  Optics  and 
Photonics  in  Global  Homeland  Security  IV,  Proceedings  of  SPIE2,  Volume  6945,  pp. 
694503  1-9,  2008. 

6.1 .2  Presentations 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Electrical  Characterization  of 
Bismuth  Iodide  Crystals,”  2010  IEEE  Nuclear  Science  Symposium  and  Medical 
Imaging  Conference,  Knoxville,  TN,  November  1-6,  2010. 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Characterization  of  Bil3  Single 
Crystals  for  Wide  Band-Gap  Semiconductor  Radiation  Detectors,”  2010  Symposium 
on  Radiation  Measurements  and  Applications  (SORMA),  Ann  Arbor,  Ml,  May  24-28, 
2010. 

•  J.  E.  Baciak,  J.  C.  Nino,  W.  Qiu,  and  A.  T.  Lintereur,  “Single  Crystal  Bismuth  Iodide 
(Bil3)  Radiation  Detectors,”  2010  ANS  Radiation  Protection  and  Shielding  Division 
Topical  Meeting,  Las  Vegas,  NV,  April  18-23,  2010. 
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•  J.  E.  Baciak,  “High  Z,  Wide  Bandgap  Semiconductor  Materials  for  Radiation 
Detection,”  2010  ANS  Student  Conference,  Ypsilanti,  Ml,  April  8-1 1 , 2010. 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Characterization  of  BH3  Single 
Crystals  Grown  by  Modified  Bridgman  Method,”  2009  IEEE  Nuclear  Science 
Symposium  and  Medical  Imaging  Conference,  Orlando,  FL,  October  25-31, 2009. 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Growth  and  Characterization  oi 
BiU  Single  Crystals  for  Room  Temperature  Gamma-Ray  Detection  Grown  by 
Modified  Bridgman  Method,”  presented  at  2009  SPIE  Optics  and  Photonics 
Conference,  San  Diego,  CA,  August  2-7,  2009. 

•  W.  Qiu,  A.  T.  Lintereur,  H.  Xu,  S.  R.  Phillpot,  J.  E.  Baciak,  and  J.  C.  Nino,  “Crystal 
Growth  and  Characterization  of  Bismuth  Tri-Iodide  by  Modified  Bridgman  Method,” 
2009  Materials  Research  Society  Spring  Meeting,  San  Francisco,  CA,  April  13-17, 
2009. 

•  W.  Qiu,  A.  T.  Lintereur,  J.  C.  Nino,  and  J.  E.  Baciak,  “Growth  and  Characterization  oi 
Single  Crystal  Bil3,”  presented  at  2009  ANS  Student  Conference,  Gainesville,  FL 
April  2-4,  2009. 

•  Azaree  T.  Lintereur,  Wei  Qiu,  Juan  C.  Nino,  and  James  E.  Baciak,  “Bismuth  Tri- 
Iodide  for  Room  Temperature  Gamma-Ray  Spectrscopy,”  2008  IEEE  Nuclear 
Science  Symposium  and  Medical  Imaging  Conference,  Dresden,  GERMANY 
October  19-25,  2008. 

•  Wei  Qiu,  Azaree  T.  Lintereur,  James  E.  Baciak,  and  Juan  C.  Nino,  "Bismuth  Iodide 
Single  Crystals  for  Gamma-Ray  Detection",  Sensors  and  Actuators  Symposium  F, 
Electroceramics  XI  conference,  Manchester,  UNITED  KINGDOM,  September  1st 
2008. 

•  W.  Qiu,  A.  T.  Lintereur,  J.  C.  Nino,  and  J.  E.  Baciak,  “Growth  and  Characterization  of 
Bismuth  Iodide  Crystals  Growth  by  the  Bridgman  Method,”  2008  SPIE  Defense  and 
Security  Conference,  Orlando,  FL  March  17-20,  2008. 

•  A.  T.  Lintereur,  W.  Qiu,  J.  C.  Nino,  and  J.  E.  Baciak,  “Iodine  Based  Compound 
Semiconductors  for  Room  Temperature  Gamma-Ray  Spectrscopy,”  2008  SPIE 
Defense  and  Security  Conference,  Orlando,  FL  March  17-20,  2008. 

•  J.  E.  Baciak,  “Iodine  Based  Compound  Semiconductor  Radiation  Detectors,”  2008 
ANS  Student  Conference,  College  Station,  TX,  February  28-March  1, 2008. 

6.2  Transitions  and  Future  Funding 

In  2007,  we  proposed  related  research  to  the  Department  of  Energy’s  NNSA  FY08  PDP 

BAA  to  investigate  electrode-surface  interfaces  for  Bil3  and  the  potential  utilization  of 

coplanar  grid  and  pixelated  electrodes.  In  April  2009,  this  proposal  was  selected  for 
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funding.  The  funding  began  on  approximately  October  1 , 2009  and  will  fund  the  Bib 
project  for  three  years,  beyond  the  DTRA  funding  ending  on  May  15,  2010. 
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